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CONVERSION  FACTORS,  NON-SI  TO  Si  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (met¬ 
ric)  units  as  follows: 


_ Multiply _ 

degrees  (angle) 

feet 

inches 

kips  (force) 

pounds  (force)  per  square  foot 
pounds  (force)  per  square  inch 
pounds  (mass)  per  cubic  foot 

pounds  (mass)  per  cubic  inch 

tons  (force)  per  square  foot 


_ ^ _ 

0.01745329 

0.3048 

25.4 

4.448222 

47.88026 

0.006894757 

16.01846 

27.6799 

95.76052 


_ To  Obtain _ 

radians 

metres 

millimetres 

kilonewtons 

pascals 

megapascals 

kilograms  per  cubic 
metre 

grams  per  cubic 
centimetre 

kilopascals 
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METHODS  OF  EVALUATING  THE  STABILITY  AND  SAFETY  OF  GRAVITY 


EARTH-RETAINING  STRUCTURES  FOUNDED  ON  ROCK 


Phase  2  Study 


PART  I :  INTRODUCTION 

1.  This  report  describes  the  continued  research  investigation  of  the 
behavior  of  earth- retaining  structures  founded  on  rock.  This  program  of  re¬ 
search  was  initiated  after  a  number  of  existing  structures,  e.g.  navigation 
lock  walls  under  the  jurisdiction  of  the  US  Army  Corps  of  Engineers,  that 
showed  no  signs  of  instability  or  substandard  performance  failed  to  meet  the 
criteria  used  for  the  design  of  new  structures.  The  traditional  analysis 
procedure  for  these  structures  uses  conventional  equilibrium  methods,  the  type 
used  widely  for  analysis  of  earth-retaining  structures.  Because  the  condi¬ 
tions  of  equilibrium  are  insufficient  for  a  complete  analysis  of  all  aspects 
of  soil -structure  interaction,  these  conventional  equilibrium  methods  neces¬ 
sarily  involve  assumptions  regarding  the  loading  and  resisting  forces  that  act 
on  the  structures.  The  finite  element  method  of  analysis  was  used  to  compare 
the  results  of  the  finite  element  analysis  with  the  conventional  equilibrium 
analysis  as  a  means  of  determining  whether  current  criteria  and  methods  are 
more  conservative  than  necessary.  The  study  of  behavior  of  earth-retaining 
structures  founded  on  rock  was  divided  into  two  phases.  The  first  phase 
(Ebeling  et  al.  in  preparation)  defined  the  analytical  requirements  for  this 
category  of  structure  and  developed  new  finite  element  procedures.  The  analy¬ 
tical  procedures,  described  in  Part  IV  of  the  first  phase  of  study,  can  be 
used  for  analysis  of  walls  that  are  loaded  so  heavily  that  gaps  develop 
between  the  base  and  the  foundation,  i.e.  structures  deemed  to  be  in  a  condi¬ 
tion  of  incipient  instability.  In  this  report,  the  second  phase  of  the  study, 
finite  element  analyses  and  equilibrium  methods  utilizing  the  results  of 
finite  element  analyses  were  performed  to  determine  the  effects  of  various 
geometrical  and  material  parameters. 

2.  In  the  finite  element  "backfill"  analyses,  the  loadings  on  the  wall 
are  generated  automatically  by  simulating  placement  of  the  backfill  behind  the 
wall.  The  magnitude  of  the  forces  acting  on  the  wall  is  dependent  upon  the 
relative  movement  of  the  soil  and  the  monolith.  These  backfill  placement 
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analyses  are  believed  to  be  the  most  realistic  that  can  be  performed  using  the 
finite  element  method. 

3.  The  finite  element  method  differs  from  the  conventional  method  of 
analysis  in  several  ways.  Two  of  the  principal  differences  are:  (a)  the 
deformations  of  the  wall  are  considered  in  the  analysis,  and  (b)  no  assump¬ 
tions  with  regard  to  the  applied  earth  loads  are  required. 

4.  Several  observations  were  made  in  the  early  studies  described  by 
Ebeling  et  al.  (in  preparation).  Among  these  were: 

a.  The  simulated  placement  of  backfill  against  the  wall  using  the 
finite  element  method  resulted  in  the  development  of  stabilizing 
shear  forces  acting  on  the  back  of  the  wall. 

b.  For  very  small  values  of  wall  displacement,  i.e.  less  than 

1  in.*  for  a  40-ft-high  wall,  the  resulting  lateral  earth  pres¬ 
sures  were  appreciably  less  than  the  at-rest  pressures. 

c.  The  lateral  translation  of  the  wall  away -from  the  backfill,  dis¬ 
placing  the  toe- fill,  resulted  in  earth  pressures  greater  than 
the  at-rest  pressures  on  the  front  of  the  wall. 

d.  As  a  result  of  these  factors,  comparisons  between  the  results  of 
backfill  placement  analyses  and  conventional  equilibrium  analy¬ 
ses  of  retaining  structures  indicated  that  the  conventional 
analysis  is  very  conservative. 

5.  The  initial  finite  element  studies  (Ebeling  et  al.  in  preparation) 
showed  that  the  shear  force  acting  on  the  back  of  the  wall  was  a  near  constant 
value  for  variations  in  the  unit  weight  and  Poisson's  ratio  of  the  backfill 
and  variations  in  the  shear  stiffness  assigned  to  the  concrete- to -rock  inter¬ 
face  along  the  base  of  the  wall.  Additional  analyses  have  been  performed  to 
assess  the  influences  on  the  shear  load  of  variations  in  these  factors: 

(a)  the  stiffness  of  the  backfill,  (b)  the  concrete-to-soil  interface  shear 
stiffness,  and  (c)  the  geometry  of  the  wall.  In  addition,  analyses  were  also 
conducted  to  study  the  effects  of  submerged  backfill  and  water  pressures  on 
the  wall. 

6.  Differences  in  the  earth  loads  calculated  in  finite  element  analyses 
(FEA)  and  those  assumed  in  conventional  equilibrium  analyses  (CEA)  reflect  the 
mobilization  of  different  levels  of  resistance  within  the  backfill  and  differ¬ 
ent  orientations  of  the  principal  stress  axes.  Methods  for  including  shear 
loads  in  equilibrium  analyses  and  tor  evaluating  the  safety  of  retaining  walls 


*  A  table  of  factors  for  converting  non- SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  5. 
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were  investigated,  as  described  in  subsequent  sections  of  this  report. 

7.  The  interdependence  between  wall  deformations  and  the  forces  acting 
on  the  wall  were  examined  in  a  series  of  equilibrium  analyses.  In  these  anal¬ 
yses,  the  deformation  of  the  wall  during  backfilling  was  sufficient  to  break 
the  cohesive  bond  between  the  wall  and  the  rock  foundation,  resulting  in  a 
situation  where  resistance  to  sliding  is  provided  only  by  friction.  Using  the 
results  of  finite  element  analyses,  relationships  between  the  height  of 
backfill,  the  earth  loads,  and  the  deformation  of  the  wall  were  developed. 

8.  This  report  is  divided  into  four  parts  and  two  appendixes.  Part  II 
discusses  the  results  of  the  advanced  backfill  placement  analyses.  These 
analyses  were  performed  to  investigate  the  effects  of  additional  variations  in 
material  parameters,  additional  wall  configurations,  and  the  presence  of  water 
behind  and  in  front  of  the  wall. 

9.  Part  III  describes  an  assessment  of  the  methods  used  for  inclusion 
of  the  factor  of  safety  in  conventional  design  analyses  of  earth-retaining 
structures . 

10.  The  findings  of  the  studies  discussed  herein  are  summarized  in 
Part  IV,  together  with  recommendations  of  topics  for  further  study. 

11.  In  Appendix  A,  two  methods  for  the  interpretation  of  the  finite 
element  results  of  the  simulated  backfill  placement  analyses  are  described. 

12.  Appendix  B  contains  a  user's  guide  for  the  modified  version  of  the 
finite  element  program  SOILSTRUCT  described  in  Part  IV  of  the  first  phase  of 
study  (Ebeling  et  al.  in  preparation). 
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PART  II:  ADVANCED  BACKFILL  PLACEMENT  ANALYSES  OF  RETAINING 
STRUCTURES  USING  THE  ALPHA  METHOD 


13.  A  series  of  advanced  backfill  placement  analyses  of  gravity- 
retaining  structures  founded  on  rock  and  performed  using  the  finite  element 
method  are  described  in  this  part.  In  a  backfill  placement  analysis,  both  the 
wall  and  the  soil  backfill  are  represented  in  the  finite  element  mesh.  The 
loadings  exerted  by  the  backfill  on  the  wall  are  generated  automatically  dur¬ 
ing  placement  of  the  backfill  behind  the  wall.  This  soil-structure  interac¬ 
tion  method  of  analysis  is  believed  to  afford  the  most  realistic  model  of  the 
backfill  procedure  for  a  wall  that  can  be  performed  using  the  finite  element 
method. 

14.  The  backfill  placement  analyses  described  in  Part  II  consider  addi¬ 
tional  variations  in  material  parameters  beyond  those  discussed  by  Ebeling  et 
al.  (in  preparation)  in  the  first  phase  of  this  study.  They  also  include 
additional  wall  configurations  and  water  behind  and  in  front  of  the  wall.  The 
finite  element  program  SOILSTRUCT,  with  the  alpha  method  incorporated  in  it, 
was  used  in  the  analyses.  The  alpha  method  allows  accurate  analysis  of  step- 
by-step  separation  of  the  base  of  t)ie  wall  from  its  foundation  during  back¬ 
filling,  as  discussed  in  the  first  phase  of  the  study.  For  those  analyses 
with  water,  uplift  pressures  were  applied  when  a  crack  developed  along  the 
base  of  the  wall. 


Structures  Analyzed 

15.  The  four  hypothetical  structures  discussed  in  this  part  are 
shown  in  Figures  1-5.  All  of  the  monoliths  retain  40  ft  of  backfill 

(H  =  40  ft),  are  16  ft  wide  at  the  base  (B  =  16  ft),  and  are  buttressed  by 
17.8  ft  of  toe -fill.  The  four  structures  differ  in  the  widths  at  the  top  of 
the  wall  and  the  slopes  along  the  back  of  the  wall.  Their  features  are  repre¬ 
sentative  of  the  existing  walls  discussed  in  Part  III  of  the  initial  study 
(Ebeling  et  al.  in  preparation). 

16.  Figure  1  shows  the  base  case  structure  used  exclusively  in  the 
backfill  placement  analyses  described  in  the  first  phase  report.  The  width  at 
the  top  of  wall  is  8  ft,  one -half  the  width  of  the  base.  No  water  table  was 
present  in  this  analysis . 
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Figure  1.  Base  case  hypothetical  structure  used 
in  the  backfill  placement  analyses 


Figure  2.  Hypothetical  structure  used  in  the  backfill 
placement  analysis  with  partially  submerged  soil 
behind  the  wall  and  a  pool  in  front 
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Very  Stiff  Rock 

Figure  3.  Rectangular  hypothetical  structure  used 
in  the  backfill  placement  analysis 


Very  Stiff  Rock 

Figure  4.  Sloping  face  hypothetical 
structure  used  in  the  backfill 
placement  analysis 


Figure  5.  Stepped  face  hypothetical  structure 
used  in  the  backfill  placement  analysis 

17.  The  same  structure  is  shown  in  Figure  2  retaining  a  partially  sub¬ 
merged  soil  backfill  and  buttressed  by  a  submerged  toe -fill  and  a  pool  of 
water.  After  completion  of  backfilling,  both  the  water  table  within  the  back¬ 
fill  and  the  pool  of  water  in  front  of  the  wall  were  at  a  height  equal  to 
26.7  ft  (0.67H).  In  subsequent  analyses,  the  elevation  of  the  pool  was 
lowered  in  two  4. 45- ft  increments  to  the  elevation  of  the  top  of  the  toe -fill. 

18.  A  rectangular  earth-retaining  wall  is  shown  in  Figure  3.  Figure  4 
shows  a  sloping  face  wall  with  a  crest  width  equal  to  1.78  ft,  or  0.11B.  A 
stepped  face  wall  with  the  same  top  width  is  shown  in  Figure  5.  There  was  no 
water  behind  or  in  front  of  the  wall  shown  in  Figures  3-5. 

Analyses  Performed 

19.  A  series  of  backfill  placement  analyses  were  performed  to  study  the 
influence  of  five  geometrical,  material,  and  loading  parameters.  The  ranges 
of  values  considered  for  the  variables  are  shown  in  Table  1.  The  purpose  of 
these  analyses  was  to  supplement  the  backfill  placement  analyses  in  the  first 
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Table  1 

Values  of  Parameters  Used  in  Backfill  Placement  Analyses 
a.  Description  of  Parameters  Used  in  Base  Case  Analysis. 


_ Constant  Parameters _ 

H  -  height  of  wall  -  40  ft 
B  -  width  of  base  of  wall  -  16  ft 
B/H  -  0.4 

htoe  "  height  of  soil  backfill  beyond 
toe 

-  17.8  ft 

htoe/H  -  0.45 

7-  “  unit  weight  of  concrete 

-  150  pcf 

Ec  “  modulus  of  elasticity  of 
concrete 
•  3  X  106  psi 

K„  *»  normal  stiffness  of  interface* 

=  3  x  106  pci 

Ks  =  shear  stiffness  of  interface* 
-lx  10'1  pci 

Poisson's  ratio  of  concrete  =0.2 

_ Variable  Parameters _ 

Ks  =  shear  stiffness  of  interface** 

=  15  to  400  pci 

Hw  =  height  of  water  behind  wall 
AH„  =  change  in  pool  elevation 
T  =  Width  at  top  of  wall 


Constant  Soil  Backfill  Parameters 

7backfiu  “  unit  weight  of  backfill 
-  135  pcf 

^backfill  “  Poisson's  ratio  of 
backfill  =0.15 

j>  =  friction  angle  of  backfill 
=  39  deg 

n  =  backfill  stiffness  parameter 
exponent  =0.4 

Rf  =  backfill  failure  ratio 
parameter  =0.7 

8  =  wall-to-soil  interface  friction 

8/ <j>  -  0.8** 


Variable  Soil  Backfill  Parameters 

K  =  backfill  stiffness  parameter 
constant 


Equations  for  Soil  Stress-Strain  Model 

Et  =  tangent  modulus  =  [1  -  R£SL]2KPa(a3/Pa)n 

SL  =  stress  level  =  ( a1  -  cr3)/(a1  -  a3)f 

(a i  -  a3)f  =  (2c  cos  <f>  +  2a3  sin  <£)/(  1  -  sin  <}>) 

P„  =  atmospheric  pressure 

(Continued) 


*  Concrete-to-rock. 

**  Concrete-to-soil. 
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Table  1  (Concluded) 

b.  Parameters  Changed  from  Values  Used  in  the  Base  Case  Analysis. 


Caset 

T 

ft 

T/B 

Hw 

ft 

< 

Backfill 

Stiffness 

Parameter 

K 

Comments ti 

17a 

8 

0.5 

-- 

450 

Base  case,  AM 

20a 

8 

0.5 

-- 

-- 

[200] 

AM 

21a 

8 

0.5 

-- 

-- 

[1,200] 

AM 

22a 

8 

0.5 

-- 

-- 

[2,500] 

AM 

23a 

8 

0.5 

-- 

-- 

[4,000] 

AM 

24a 

8 

0.5 

-- 

-- 

[6,000] 

AM 

25a 

8 

0.5 

-- 

-- 

450 

[(Ks**)  1,  AM 

min 

26a 

8 

0.5 

-- 

450 

[(Ks**)  ],  AM 

max 

27b 

8 

0.5 

[26.7] 

[0] 

450 

AM 

28b 

8 

0.5 

[26.7] 

[4.45] 

450 

AM 

29b 

8 

0.5 

[26.7] 

[8.9] 

450 

AM 

30a 

[16] 

1.0 

-- 

450 

AM 

31a 

[1.8] 

0.11 

-- 

-- 

450 

AM 

32a 

[1.8] 

0.11 

450 

[stepped  wall] ,  . 

Note:  --  indicates  no  data.  Parameters  in  brackets  are  those  changed  from 

the  values  for  the  base  case. 

**  Concrete-to-soil. 

f  a  -  no  hydrostatic  uplift  pressure  applied  to  area  of  separation  along 
the  base,  b  =  full  hydrostatic  uplift  pressure  applied  to  area  of 
separation  along  the  base, 
ft  AM  =  alpha  method. 
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phase  of  study  (Ebeling  et  al.  in  preparation)  by  investigating  the  influence 
of  the  following  parameters  on  the  calculated  results: 

a.  The  stiffness  assigned  to  the  backfill. 

b.  The  interface  stiffnesses. 

c.  The  presence  of  water  in  front  of  and  behind  the  wall. 

d.  Variations  in  wall  geometry. 

20.  As  in  previous  backfill  placement  analyses,  a  number  of  factors 
that  affect  the  stability  of  the  wall  were  studied.  These  include: 

a.  The  distribution  of  stresses  along  the  base  of  the  wall,  the 
front  and  back  of  the  wall,  and  along  several  planes  through 
the  backfill. 

b.  The  magnitudes  of  the  resultant  forces  on  these  planes,  and  the 
positions  of  their  points  of  action. 

c.  The  magnitudes  of  the  mobilized  angles  of  friction  on  these 
planes . 

d.  The  magnitudes  of  the  earth  pressure  coefficients  that  char¬ 
acterize  the  magnitudes  of  the  earth  pressures  on  vertical 
planes  within  the  backfill. 

e.  The  percent  of  effective  contact  between  the  base  of  the  wall 
and  the  foundation  (Be/B) . 

f.  The  lateral  displacement  of  the  monolith. 

g.  The  magnitude  of  the  maximum  compressive  stress  developed  at 
the  toe  of  the  wall. 

Effect  of  the  Stiffness  of  the  Backfill 

21.  Six  analyses  were  performed  to  investigate  the  influence  of  back¬ 
fill  stiffness.  The  base  case  structure  shown  in  Figure  1  was  used.  The 
stiffness  of  the  backfill  was  varied  to  cover  a  range  of  values  characteristic 
of  typical  backfill  soils. 

Loading  scheme 

22.  In  backfill  placement  analyses,  the  loadings  exerted  on  the  wall  by 
the  backfill  are  dependent  upon  the  magnitude  and  direction  of  wall  movement. 
These  loads  are  not  predetermined,  but  depend  on  soil-structure  interaction. 

23.  Prior  to  backfilling,  the  construction  of  the  wall  was  simulated  in 
three  lifts,  each  13.33  ft  high.  The  simulation  is  idealized  in  Figure  6  for 
the  base  case  structure.  Construction  of  the  wall  was  followed  by  placement 
of  40  ft  of  backfill  in  nine  layers,  each  4.44  ft  thick.  The  first  four 
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Concrete  Bockfill 

Lift  Layer 


Figure  6.  Idealization  of  the  base  case 
structure  for  backfill  placement 
analysis- -no  water 


layers  were  placed  simultaneously  in  front  of  and  behind  the  wall  to  a  height 
of  17.8  ft.  No  water  pressures  were  represented  in  this  analysis. 

Finite  element  mesh 

24.  The  finite  element  mesh  used  for  the  backfill  placement  analysis  of 

the  base  case  structure  is  shown  in  Figure  7.  This  mesh  and  the  other  meshes 
used  in  these  analyses  were  designed  to  model  gravity  walls  founded  on  very 
stiff  rock.  The  nodes  along  the  base  were  fixed,  and  thus  simulated  a  rigid 
boundary  at  the  top  of  rock.  Interface  elements  were  included  in  four 
regions:  (a)  between  the  wall  and  the  backfill,  (b)  between  the  wall  and  the 

foundation,  (c)  between  the  backfill  and  the  rock,  and  (d)  between  the  toe- 
fill  and  the  rock. 

25.  There  are  258  two-dimensional  and  interface  elements  in  the  mesh 
shown  in  Figure  7;  57  elements  model  the  wall,  117  elements  model  the  back¬ 
fill,  and  40  elements  model  the  toe -fill.  The  remaining  elements  model  the 
interfaces . 
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Figure  7 .  Finite  element  mesh  used  to  model  a  base  case  hypothetical 
structure  with  additional  backfill  beyond  the  toe 


Material  properties 

26.  As  in  the  backfill  placement  analyses  discussed  in  Pa^t  VI  of  the 
first  phase  of  this  study  (Ebeling  et  al.  in  preparation),  the  wall  was 
modeled  as  a  linear  el^jtic  material  and  the  soil  as  a  hyperbolic,  nonlinear, 
stress -dependent  material.  The  values  of  the  parameters  used  in  these  analy¬ 
ses  are  listed  in  Table  1  along  with  some  of  the  key  equations  of  the  hyper¬ 
bolic  stress-strain  relationships. 

27 .  The  features  of  the  interfaces  were  discussed  in  Part  VI  of  the 
initial  phase  of  this  study  by  Ebeling  et  al.  (in  preparation).  A  very  useful 
feature  of  the  interface  along  the  base  of  the  wall  is  it^  ability  to  model 
the  development  of  a  crack  during  backfilling,  using  the  alpha  method  (AM)  of 
analysis  described  in  Part  IV  of  the  Ebeling  et  al.  (in  preparation)  study. 

28.  The  properties  of  the  concrete  and  the  interface  between  the  wall 

and  its  foundation  used  in  these  analyses  are  the  same  as  those  used  in  the 

backfill  placement  analyses  discussed  in  the  first  phase  of  study.  The  modu¬ 
lus  of  the  concrete  was  3,000,000  psi,  and  the  Poisson's  ratio  of  the  concrete 
was  0.15.  The  normal  stiffness  of  the  interface  between  the  wall  and  the 
foundation  was  3,000,000  pci,  and  the  shear  stiffness  was  10,000  pci. 

29.  In  the  backfill  analyses  discussed  initially  (Ebeling  et  al.  in 

preparation) ,  the  properties  assigned  to  the  backfill  and  the  interfaces  were 
characteristic  of  a  clean  granular  backfill  with  a  relative  density  of  about 
75  percent.  The  unit  weight  was  135  pcf,  and  the  angle  of  internal  friction 

was  39  deg.  The  magnitudes  of  the  parameters  used  to  model  the  hyperbolic 

stress-strain  behavior  were:  the  modulus  number,  K  **  450  ;  the  modulus  expo¬ 
nent,  n  «=  0.4  ;  and  the  failure  ratio,  Rf  =  0.7  .  The  typical  range  in 
values  of  K  for  clean  granular  backfill  is  from  200  to  1,200.  In  this 
series  of  six  analyses,  the  value  of  K  was  varied  from  a  minimum  value  of 

K  -  200  to  a  maximum  value  of  K  =  6,000  ,  with  intermediate  values  of 
K  -  450  (base  case  analysis),  1,200,  2,500,  and  4,000.  Values  of  K  greater 
than  1,200  were  used  in  the  analyses  so  that  the  trends  in  the  results  for 
extreme  values  of  K  could  be  established. 

30.  The  shear  stiffness  values  assigned  to  the  interfaces  between  the 
soil  and  the  concrete  and  between  the  soil  and  the  rock  ranged  from  20  to 
300  pci,  depending  on  the  confining  pressure.  The  normal  stiffness  assigned 
to  all  of  the  interfaces  was  3,000,000  pci.  The  friction  angle  of  the 
wall/soil  interface  was  31  deg,  about  80  percent  of  the  angle  of  internal 
friction  for  the  soil. 
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Results  of  the 

backfill  placement  analyses 

31.  The  results  of  the  analyses  (Cases  17a,  20a,  21a,  22a,  23a,  and 
24a)  are  summarized  in  Table  2.  The  resultants  of  the  stress  distributions 
computed  using  the  program  SOILSTRUCT  were  computed  manually.  Their  magni¬ 
tudes  and  points  of  action  were  determined  for  each  of  the  several  planes 
shown  in  Figure  8.  The  Sections  labeled  A-A  and  C-C  pass  through  the  heel  of 
the  monolith.  Section  A-A  extends  vertically  within  the  backfill,  and  Sec¬ 
tion  C-C  passes  along  the  soil-to-wall  interface.  Section  E-E  is  a  vertical 
plane  at  the  toe  of  the  monolith.  Vertical  Sections  B-B  and  H-H  are  located 
far  behind  and  far  in  front  of  the  wall.  The  resultant  forces  acting  on  these 
planes  were  used  to  determine  the  at-rest  earth  pressure  coefficients,  which 
characterize  the  lateral  loads  on  an  unyielding  wall. 

Forces  on  Section  B-B 

32.  Due  to  its  distance  from  the  monolith,  the  stresses  which  developed 
on  the  vertical  Section  B-B  were  not  influenced  by  soil-structure  interaction. 
The  stresses  on  this  plane,  therefore,  reflect  only  the  influence  of  the  mate¬ 
rial  parameters  assigned  to  the  soil. 

33.  The  magnitude  of  the  resultant  lateral  force  acting  on  Section  B-B 
F^  was  found  to  be  nearly  exactly  the  same  for  all  values  of  the  modulus 
number  K  .  For  example,  with  K  equal  to  450  (base  case),  F^  was  computed 
to  be  55,329  lb,  and  with  K  equal  to  the  maximum  value  of  6,000,  F^  was 
found  to  be  55,137  lb,  a  difference  of  only  192  lb,  or  less  than  1  percent. 

34.  The  lateral  earth  pressure  coefficient  for  Section  B-B  is 

equal  to  the  at-rest  earth  pressure  coefficient  KQ  .  Values  of  ■*  Kq 
were  calculated  from  the  finite  element  results  using  the  expression 


K.  „  _ — _ 

“  /"H  -  « 

/  a  dh 

Jo  v 

where 

F^  °  resultant  lateral  effective  force 
H  -  height  of  soil  above  the  base  of  the  wall 
av  “  the  vertical  effective  stress  at  height  h 


(1) 
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Table  2 

Summary  of  Results  for  Backfill  Placement  Analyses  of  the  Base  Case  Structure 
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(Concluded) 
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Concrete 

Monolith 


Planes  along  surfaces  of  the  monolith  and  within  the  backfill 
on  which  resultant  forces  are  computed- -no  water 


With  no  water  in  the  backfill,  Equation  1  simplifies  to  the  relationship 


^backfill  ‘  H 


(2) 


For  the  six  analyses  summarized  in  Table  2,  KQ  is  equal  to  0.51  using  Equa¬ 
tion  2.  The  value  of  KQ  serves  as  a  convenient  index  to  the  magnitude  of 
loading  this  wall  would  be  subjected  to  if  it  did  not  move. 

35.  It  may  be  seen  that  the  resultant  lateral  force  acts  at  0.36H  above 
the  base  in  all  six  analyses.  This  is  slightly  higher  than  the  value  of  0.33H 
corresponding  to  an  exactly  triangular  at-rest  pressure  distribution. 

36.  The  values  of  Kv  shown  in  Table  2  were  calculated  using  the 
equation: 


K 


v  r  h 
0 


a  dh 
v 


(3) 


The  resulting  vertical  shear  force  on  Section  B-B  (F  )  is  zero  for  all  analy¬ 
ses  as  a  result  of  the  fact  that  there  are  no  differential  settlements  from 
one  side  of  Section  B-B  to  the  other.  The  vertical  shear  stress  coefficient 
(Kv)  at  this  section  is,  thus,  zero,  since  Fv  is  zero.  For  the  analyses 
where  no  water  table  is  present,  Equation  3  becomes 


K 

v 


0.5 


F 

v 


"^backfill 


(4) 


37.  The  values  of  Sm  shown  in  Table  2  were  calculated  using  the 
equation: 


tan  (5  ) 
m 


F 


v 


The  mobilized  angle  of  friction  (5  )  is  also  zero  at  Section  B-B. 


(5) 
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Forces  on  the  wall 


38.  The  values  of  the  resultant  forces  on  the  monolith  are  summarized 
in  Table  2.  The  magnitude  of  the  shear  force  acting  on  the  base  of  the  mono¬ 
lith  (T)  decreased  with  increasing  values  of  the  modulus  number  K  .  However, 
for  values  of  K  within  the  range  for  typical  backfill  soils  and  values  of  K 
between  200  and  1,200,  T  was  nearly  constant  with  a  value  of  37,400  lb. 

39.  The  resultant  normal  force  on  the  base  of  the  wall  (N)  decreased 
slightly  as  the  value  of  K  increased.  An  increase  in  the  value  of  K  from 
200  to  1,200  resulted  in  a  4.5  percent  decrease  in  the  value  of  H  .  It  may 
be  seen  that  distance  between  the  resultant  force  and  the  toe  of  the  wall  (x  ) 
decreased  as  the  value  of  K  increased.  This  reflects  the  decrease  in  the 
magnitude  of  the  shear  i  res s  applied  to  the  wall  by  the  backfill  as  the 
stiffness  of  the  backfill  increased. 

40.  The  values  of  x*  shown  in  Table  2  are  the  distances  from  the  toe 
to  the  resultants  of  the  normal  force  for  the  elements  remaining  in  compres¬ 
sion.  It  may  be  seen  that  there  is  good  agreement  between  the  values  of  xn 
and  x*  ,  indicating  that  the  numerical  accuracies  of  the  results  are  very 
good. 

41.  The  magnitude  of  the  normalized  effective  base  contact  area  (B  /B) 
decreased  with  increasing  backfill  stiffness.  Fifty-six  percent  of  the  base 
remained  in  compression  for  a  value  of  K  equal  to  200,  while  the  effective 
base  contact  area  was  reduced  to  50  percent  for  a  value  of  K  equal  to  1,200. 

42.  As  the  stiffness  of  the  backfill  increased,  the  maximum  compressive 
stress  at  the  toe  of  the  monolith  (q__„)  increased.  As  the  value  of  K  was 
increased  within  the  range  of  values  representative  of  backfill  soils,  9max 
increased  from  32,106  to  34,367  psf,  an  increase  of  less  than  1  percent. 

43.  The  mobilized  angle  of  internal  friction  along  the  base  (5^) 
increased  with  increasing  backfill  stiffness.  As  the  backfill  modulus  number 
K  increased  from  200  to  6,000,  the  value  of  6'm^  increased  from  18.5  to 
20.4  deg.  The  change  results  from  the  decrease  in  the  magnitude  of  the  normal 
force  on  the  base,  since  the  shear  force  is  nearly  constant. 

44.  The  lateral  movement  of  the  monolith  also  increased  with  increasing 
values  of  K  .  The  normalized  lateral  deformation  at  the  crest  of  the  mono¬ 
lith  (ux/H)  increased  from  0.000049  for  K  =  200  to  0.000080  for  K  =  6,000  . 
These  values  of  wall  deformation  are  extremely  small  compared  to  the  values 
which  Terzaghi  (1934)  found  were  needed  for  the  development  of  active  pres¬ 
sures  in  his  full  scale  tests  on  dense  sand  (u^/H  »  0.0014). 
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Forces  on  Section  C-C 


45.  Section  C-C  is  the  interface  between  the  wall  and  the  backfill. 

The  forces  acting  on  this  plane  are  affected  by  the  value  of  KQ  for  the 
backfill,  the  settlement  of  the  backfill  as  it  is  placed,  and  the  effects  of 
the  movements  of  the  wall  during  backfilling. 

46.  The  magnitude  of  the  resultant  shear  force  (F  )  decreased  with 
increasing  stiffness  of  the  backfill.  For  K  =  200  ,  the  value  of  Fs  was 
25,833  lb,  and  for  K  -  1,200  ,  the  value  of  Fg  was  22,739  lb,  a  decrease  of 
3,094  lb  (12.0  percent).  This  is  due  to  a  reduction  in  the  amount  of  settle¬ 
ment  of  the  backfill  as  the  soil  stiffness  increased. 

47.  The  normal  force  on  Section  C-C  (F  )  remained  nearly  constant, 
varying  only  from  54,681  to  55,906  lb  as  K  was  varied.  For  all  six  analy¬ 
ses,  the  resultant  normal  force  acted  at  a  distance  of  14  ft  above  the  heel, 
as  measured  along  the  face.  This  length  corresponds  to  35  percent  of  the 
interface  length  (L) . 

48.  The  mobilized  angle  of  friction  for  the  interface  (6  )  decreased 
from  25.3  deg  for  K  =  200  to  24.5  deg  for  K  =  450  ,  and  to  22.2  deg  for 

K  *=  1,200  .  The  decrease  in  Sm  with  increasing  values  of  K  results  from  a 
decrease  in  the  settlement  of  the  backfill  and  the  magnitude  of  F  . 

Forces  on  Section  A-A 

49.  Section  A-A  is  a  vertical  plane,  within  the  backfill,  passing 
through  the  heel  of  the  wall.  Like  Section  C-C,  the  forces  acting  on  this 
plane  are  affected  by  the  soil-structure  interaction. 

50.  The  magnitude  of  the  vertical  shear  force  (F  )  acting  on  Sec¬ 
tion  A-A  decreased  with  increasing  soil  stiffness.  Fv  was  equal  to  15,435  lb 
for  K  =  200  ,  decreasing  to  14,795  lb  for  the  base  case  value  of  K  =  450  , 
and  to  12,718  lb  for  K  -  1,200  .  Figure  9  shows  this  trend.  The  value  of 

Kv  would  be  expected  to  approach  zero  as  the  stiffness  of  the  backfill 
approaches  the  stiffness  of  concrete.  For  the  typical  range  in  stiffnesses  of 
good  quality  backfill  materials,  Kv  ranges  in  value  from  about  0.12  to  about 
0.15.  The  downdrag  force,  Fv  ,  provides  a  stabilizing  influence  on  the  wall. 

53 .  The  magnitude  of  the  resultant  lateral  force  on  Section  A-A  (F^) 
increased  with  increasing  values  of  K  .  The  lateral  earth  pressure  coeffi¬ 
cient  (K^)  was  equal  to  0.45  for  K  =  200  ,  increasing  to  0.47  for  K  ■=  6,000  . 
For  the  typical  range  in  K  values  for  good  quality  backfill,  the  value  of 
the  lateral  earth  pressure  coefficient  for  Section  A-A  (K^)  divided  by  KQ 
range  from  0.88  to  0.91. 
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Backfill  Stiffness  Parameter  K 

Figure  9.  Variation  of  vertical  earth  pressure  coefficient 
with  backfill  stiffness  parameter 


52.  The  earth  pressure  index  (EPI)  reflects  the  magnitude  of  for 

Section  A-A,  compared  to  KQ  and  Ka  .  The  EPI  is  given  by 


EPI  = 


(6) 


A  value  of  EPI  =  0  represents  a  condition  of  no  wall  movement  and  no  soil- 
structure  interaction.  Complete  soil -structure  interaction  occurs  when  a 
value  of  EPI  =1.0  is  reached.  This  represents  a  condition  where  the  wall 
movements  have  reached  their  maximum  effect,  and  the  minimum  magnitude  force 
is  exerted  by  the  backfill  on  the  wall  (K^  =  Ka) .  For  the  value  of  K 
between  200  and  1,200,  the  EPI  was  found  to  vary  from  0.17  to  0.22. 

53.  The  mobilized  angle  of  internal  friction  for  Section  A-A  ($_)  was 
found  to  decrease  with  increasing  values  of  backfill  stiffness.  It  may  be 
seen  that  Sm  decreased  from  17.6  deg  for  K  =  200  ,  to  16.8  deg  for 
K  -  450  ,  and  to  14.2  deg  for  K  =  1,200  .  This  variation  reflects  the 
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simultaneous  decrease  in  F  and  increase  in  Fu  ,  which  occur  with  increas- 

v  n 

ing  values  of  backfill  stiffness. 

Forces  on  Section  H-H 

54.  Like  Section  B-B,  the  stresses  developed  on  Section  H-H  are  not 
affected  by  soil-structure  interaction  and  reflect  only  the  parameters  of  the 
fill.  The  lateral  earth  pressure  coefficient  (K^)  at  this  section  corresponds 
to  at-rest  pressure  and  is  equal  to  0.6  for  all  six  analyses.  Since  the 
height  and  the  number  of  lifts  used  to  simulate  the  toe-fill  are  not  the  same 
as  those  used  for  the  backfill,  the  lateral  earth  pressure  coefficients  for 
these  two  sections  are  not  equal.  This  is  due  to  the  incremental  computa¬ 
tional  procedure  used  in  the  analyses,  as  discussed  in  Part  VI  of  the  first 
phase  of  this  study  by  Ebeling  at  al.  (in  preparation). 

Forces  on  Section  E-E 

55.  Section  E-E  is  a  vertical  plane  passing  through  the  toe  of  the 

monolith.  It  was  observed  that  as  the  stiffness  of  the  soil  increased,  the 
magnitude  of  the  vertical  earth  pressure  coefficient  (Kv)  decreased  and  the 
magnitude  of  the  lateral  earth  pressure  coefficient  (K^)  increased.  The  value 
of  Kv  decreased  from  0.20  for  K  =  200  to  a  value  of  0.09  for  K  “  1,200  , 
while  increased  from  0.52  to  0.60.  The  change  in  the  resultant  forces  on 

Section  E-E  are  the  combined  results  of  the  changes  in  the  magnitude  of  the 
settlement  of  the  toe-fill,  and  movement  of  the  wall  toward  the  toe-fill. 

Conclusions  on  the  effect 
of  stiffness  of  the  backfill 

56.  In  the  hyperbolic  soil  model,  two  parameters  control  the  stress - 
strain  behavior  of  the  backfill:  (a)  modulus  number,  and  (b)  Poisson's  ratio. 
It  was  observed  in  the  parametric  studies  described  by  Ebeling  et  al.  (in 
preparation)  that  variations  in  the  value  of  Poisson's  ratio  assigned  to  the 
backfill  resulted  in  variations  in  the  magnitude  of  the  lateral  force  acting 
on  the  wall.  In  this  series  of  parametric  studies,  changes  in  the  stiffness 
assigned  to  the  backfill  influenced  the  magnitude  of  both  the  shear  and 
lateral  forces  acting  on  the  wall.  As  the  stiffness  of  the  backfill 
increased,  the  magnitude  of  the  shear  stress  acting  on  the  back  of  the  wall 
decreased  slightly,  while  the  magnitude  of  lateral  force  exerted  on  the  wall 
increased  slightly.  The  downdrag  force  exerted  by  the  backfill  provides  a 
stabilizing  influence  on  the  wall,  and  a  decrease  in  its  magnitude  results  in 
a  lower  margin  of  safety  against  overturning  about  the  toe  of  the  wall. 
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Effect  of  the  Interface  Shear  Stiffness 


57 .  Two  analyses  were  performed  to  determine  the  influence  of  changing 
the  shear  stiffness  (Kg)  of  the  interfaces  between  the  concrete  and  soil  and 
between  the  rock  and  soil.  The  mesh  shown  in  Figure  7  was  used  in  these  anal¬ 
yses.  In  Case  25a,  the  values  of  Kg  used  in  the  analyses  ranged  from  17  to 
45  pci,  depending  on  the  normal  stress  in  the  elements.  In  Case  26a,  the 
values  of  Kg  ranged  from  23  to  405  pci.  The  values  of  all  of  the  other 
parameters  were  the  same  as  in  the  base  case  analysis  (17a) . 

Forces  on  Section  B-B 

58.  The  magnitude  of  the  resultant  lateral  force  (F^)  and  shear  force 

(F  )  acting  on  Section  B-B  was  the  same  for  both  analyses,  as  seen  in  Table  2. 
This  is  in  agreement  with  the  observation  that  the  magnitude  of  the  stresses 
developed  within  this  region  of  the  mesh  are  dependent  solely  on  the  proper¬ 
ties  assigned  to  the  backfill.  As  in  the  base  case  analysis,  was  equal 

to  0.51,  and  Ky  and  5m  were  equal  to  zero. 

Forces  on  the  wall 

59.  The  magnitude  of  the  resultant  shear  force  acting  along  the  base  of 
the  wall  (T)  decreased  from  39,993  to  37,121  psf  with  increasing  concrete-to- 
soil  and  rock-to-soil  interface  stiffnesses.  The  reduction  in  the  shear  force 
on  the  base  of  the  wall  is  due  to  the  fact  that  when  the  interface  between  the 
soil  and  the  underlying  rock  foundation  is  stiffer,  movement  of  the  backfill 
toward  the  wall  is  inhibited,  and  the  earth  load  on  the  wall  is  reduced. 

60.  When  Kg  was  increased,  the  magnitude  of  the  normal  force  on  the 
base  of  the  wall  (N)  increased  and  its  point  of  action  (x  )  moved  nearer  to 
the  heel  of  the  wall.  This  reflects  an  increase  in  the  magnitude  of  downdrag 
force  and  a  decrease  in  lateral  force  acting  on  the  back  of  the  wall. 

61.  The  normalized  effective  base  contact  area  (Be/B)  increased  from 
0.44  for  the  minimum  values  of  Kg  to  0.50  for  the  standard  values,  and  ulti¬ 
mately  to  0.56  for  the  maximum  values  of  Kg  used  in  the  analyses.  As  the 
effective  base  contact  area  increased,  the  maximum  compressive  stress  at  the 
toe  of  the  wall  decreased  from  37,978  to  31,711  psf. 

62.  The  mobilized  angle  of  friction  along  the  base  (6  decreased  from 
20.8  to  18.5  deg  as  the  interface  stiffnesses  increased. 

63.  The  normalized  lateral  deformation  at  the  crest  of  the  monolith 
(u/H)  decreased  from  0.000063  to  0.000049  as  a  result  of  the  decrease  in  the 
magnitude  of  the  lateral  force  applied  to  the  wall. 
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Forces  on  Section  C-C 

64.  An  increase  in  the  shear  stiffness  of  the  interface  between  the 
back  of  the  wall  and  the  backfill  resulted  in  an  increase  in  the  magnitude  of 
the  shear  force  acting  along  Section  C-C.  The  normal  force  decreased  slightly 
as  a  result  of  the  increase  in  stiffness  of  the  interface  between  the  backfill 
and  the  underlying  rock.  The  elevation  of  the  resultant  normal  force  remained 
at  the  same  point  along  the  interface  in  the  three  analyses.  The  mobilized 
friction  angle  increased  from  19.4  to  24.8  deg  as  the  interface  stiffnesses 
were  increased. 

Forces  on  Section  A-A 

65.  The  shear  force  on  vertical  Section  A-A  increased  and  the  lateral 

force  decreased  as  the  magnitude  of  Kg  increased.  The  corresponding  values 
of  Ky  ranged  in  value  from  0.09  to  0.14  and  the  values  of  from  0.48  to 

0.45.  The  points  of  action  of  the  resultant  normal  forces  in  the  three  analy¬ 
ses  differ  by  less  than  1  ft.  The  value  of  6  increased  from  11.1  to 

17.1  deg  as  Ks  was  increased.  The  value  of  the  EPI  increased  from  0.10  to 
0.21,  indicating  a  greater  level  of  soil -structure  interaction  with  larger  Kg 
values . 

Forces  on  Section  H-H 

66.  As  was  the  case  for  Section  B-B,  the  resultant  forces  acting  on 
vertical  Section  H-H  are  the  same  for  all  three  analyses.  This  is  due  to  the 
lack  of  soil-structure  interaction  at  large  distances  from  the  wall.  The 
lateral  earth  pressure  coefficient  (K^)  at  this  section  corresponds  to  Kq 
and  is  equal  to  0.6.  With  Fv  equal  to  zero,  Kv  and  Sm  are  also  equal  to 
zero  for  Section  H-H. 


Forces  on  Section  E-E 

67.  The  resultant  shear  force  on  vertical  Section  E-E  is  the  net  effect 
of  two  different  modes  of  displacement,  the  initial  settlement  of  the  toe -fill 
during  backfilling  followed  by  the  upward  movement  of  the  fill  near  the  toe  as 
the  wall  displaces  toward  the  toe -fill.  An  increase  in  the  shear  stiffness 
results  in  an  increase  in  the  magnitude  of  the  shear  forces  resulting  from 
both  of  these  displacements.  It  was  observed  that  the  wall  movements  were 


quite  small,  and  thus  the  settlement  of  the  toe- fill  during  backfilling  con¬ 


trolled  the  development  of  the  resulting  interface  forces. 


68.  The  values  of  the  lateral  force  were  nearly  constant  for  the  three 


analyses.  Expressed  in  terms  of  ,  their  values  varied  from  0.55  to  0.57 
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for  the  three  analyses.  The  interface  shear  stresses,  expressed  in  terms  of 
,  increased  from  0.12  to  0.16  as  the  interface  shear  stiffnesses  were 
increased. 

Conclusions  on  the  effect 
of  interface  shear  stiffness 

69.  Increases  in  the  values  of  the  shear  stiffness  of  the  concrete- to- 
soil  interfaces  result  in  a  more  stable  wall  because  it  results  in  larger 
shear  stresses  along  the  front  and  back  faces  of  the  wall.  Also,  changes  in 
the  value  of  Kg  between  the  backfill  and  the  underlying  rock  result  in 
reduced  earth  loads  on  the  wall,  and  thus  a  more  stable  wall. 

Effect  of  Water  in  Front  of  and  Behind  the  Wall 


70.  Three  analyses  were  performed  to  determine  the  influence  of  a  par¬ 
tially  submerged  backfill,  submerged  toe -fill,  and  a  pool  of  water  in  front  of 
the  wall.  The  base  case  retaining  structure,  shown  in  Figure  2,  was  used  in 
the  evaluation.  The  results  for  Case  27b  represent  the  completion  of  the 
backfill  placement  analysis  with  the  elevation  of  the  water  outside  the  wall 
the  same  as  that  of  the  water  table  within  the  backfill  (26.7  ft).  The  pool 
of  water  in  front  of  the  wall  was  subsequently  lowered  in  two  4.45-ft  incre¬ 
ments  to  the  elevation  of  the  top  of  the  toe-fill,  17.8  ft  above  the  base  of 
the  wall  (Cases  28b  and  29b) .  The  analyses  were  performed  using  the  finite 
element  mesh  for  the  base  case  structure  shown  in  Figure  7.  All  the  parame¬ 
ters  were  the  same  as  the  base  case  (17a)  except  that  the  buoyant  unit  weight 
was  assigned  to  the  submerged  backfill. 

Loading  scheme 

71.  Figure  10  shows  the  method  used  for  incorporating  water  pressures 
in  the  backfill  placement  analyses.  For  this  analysis  the  water  table  was  as¬ 
sumed  to  rise  as  the  backfill  was  placed.  The  simultaneous  placement  of  back¬ 
fill  layer  6  and  the  rise  in  water  table  from  the  top  of  layer  5  to  the  top  of 
layer  6  is  shown  in  Figure  10.  This  procedure  models  the  behavior  of  a  free- 
dreining  backfill.  By  assigning  buoyant  unit  weights  to  submerged  backfill 
layer  6  during  the  analysis,  the  computed  effective  stresses  within  the  back¬ 
fill  were  consistent  with  fully  drained  conditions  in  the  backfill. 

72.  Water  pressures  were  applied  normal  to  the  front  and  back  of  the 
wall  simulating  the  simultaneous  rise  in  water  table  on  either  side  of  the 
wall.  The  incremental  change  in  hydrostatic  pressures  shown  in  Figure  10 


31 


correspond  to  the  change  in  water  table  from  the  top  of  layer  5  to  the  top  of 
layer  6.  In  this  manner,  the  full  hydrostatic  water  pressures  acting  on  the 
wall  were  applied  incrementally  during  all  stages  of  backfilling. 

73.  Also  depicted  in  Figure  10  is  the  possible  separation  of  a  portion 
of  the  base  of  the  wall  from  its  foundation.  Due  to  the  development  of  such  a 
crack  along  the  base,  water  would  fill  the  void  and  exert  a  pressure  equal  to 
the  hydrostatic  pressure  in  the  backfill.  This  water  pressure  acts  on  both 
the  wall  and  the  rock  foundation  within  the  cracked  region  of  the  interface. 

74.  The  method  used  to  model  the  lowering  of  the  pool  in  front  of  the 
wall  is  shown  at  the  right  side  of  Figure  10.  For  each  of  the  two  changes  in 
the  pool  elevation  (AH„  -  4.45  ft),  a  pressure  distribution  equal  to  the 
change  in  hydrostatic  pressure  was  applied  normal  to  the  face,  directed  away 
from  the  front  of  the  wall. 

Results  of  the 

backfill  placement  analyses 

75.  The  results  of  the  three  analyses  (27b,  28b,  and  29b)  are  sum¬ 
marized  in  Table  3.  The  difference  in  the  three  cases  is  the  value  of  hu  , 
the  elevation  of  this  resultant  water  pressure  force  in  front  of  the  wall. 

The  values  of  hu  are  8.9  ft  for  Case  27b,  7.4  ft  for  Case  28b,  and  5.9  ft 
for  Case  29b.  The  analyses  were  summarized  using  resultant  effective  forces 
and  the  forces  due  to  water  pressures  acting  on  the  planes  shown  in  Figure  11, 
the  same  planes  used  in  previous  discussions. 

Forces  on  Section  B-B 

76.  With  the  same  material  parameters  and  the  same  water  table  eleva¬ 

tion  in  the  backfill  in  the  three  analyses ,  the  magnitude  of  the  resultant 
lateral  effective  force  (F^)  acting  on  Section  B-B  was  nearly  a  constant 
value  of  44,100  lb/ft.  The  resultant  water  pressure  force  (U)  was  a  constant 
22,242  lb/ft.  The  point  of  action  of  and  U  were  at  14.1  and  8.9  ft 

above  the  base  of  the  wall.  The  value  of  was  equal  to  0.51,  in  accor¬ 

dance  with  Equation  1.  With  no  differential  settlements  at  Section  B-B,  the 
resultant  shear  force  (F  )  was  equal  to  zero,  as  were  (Equation  3)  and 

5m  • 

Forces  on  the  wall 

77.  The  toe -fill  and  the  pool  of  water  in  front  of  the  wall  both  have  a 
stabilizing  influence,  because  the  base  is  not  required  to  provide  the  entire 
resistance  to  the  forces  acting  on  the  back  of  the  wall.  As  shown  in  Fig¬ 
ure  12,  the  magnitude  of  the  shear  force  along  the  base  (T)  increased  as  the 
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Forces  acting  on  the  wall  versus  change  in  pool  level  in  front 

of  the  wall 


pool  elevation  was  lowered.  Initially,  T  provided  53  percent  of  the  total 
lateral  resistance.  The  contribution  of  T  increased  to  71  percent  when  the 
elevation  of  the  pool  was  lowered  8.9  ft. 

78.  The  magnitude  of  the  effective  normal  force  (N')  decreased  with 
decreasing  elevation  of  the  water  in  front  of  the  wall,  as  a  result  of  a 
decrease  in  effective  base  contact  area  and  corresponding  increase  in  uplift 
pressures.  When  the  pool  was  lowered  8.9  ft,  the  value  of  N'  decreased  from 
89,417  to  83,297  lb,  a  decrease  of  6,120  lb.  At  the  same  time,  the  uplift 
pressure  increased  by  6,664  lb.  The  vertical  shear  forces  acting  on  the  front 
and  back  of  the  wall  accounted  for  the  difference  of  544  lb. 

79.  When  the  water  level  in  front  of  the  wall  was  lowered  by  8.9  ft, 
the  normalized  effective  base  contact  area  (Bg/B)  decreased  by  62  percent  from 
a  value  equal  to  0.34  to  a  value  of  0.13.  The  mobilized  angle  of  internal 
friction  along  the  base  (6  increased  from  19.8  to  26.9  deg.  The  normalized 
lateral  deformation  of  the  crest  of  the  wall  (ux/H)  increased  by  a  factor  of 
two  and  one -half  when  the  pool  was  lowered,  from  0.000064  to  0.000161. 

Forces  on  Section  C-C 

80.  An  8.9-ft  decrease  in  pool  elevation  in  front  of  the  wall  resulted 
in  a  5-percent  increase  in  the  magnitude  of  the  shear  force  (Fg)  and  a 

3  percent  decrease  in  the  magnitude  of  the  normal  effective  force  (F^)  acting 
on  the  back  of  the  wall.  The  change  in  the  magnitudes  of  Fg  and  F^  were 
attriDuted  to  increased  lateral  deformations  of  the  wall  away  from  the  back¬ 
fill.  The  effective  normal  force  acts  at  the  same  location  along  the  inter¬ 
face  in  all  three  cases.  The  mobilized  friction  angle  for  Section  C-C 
increased  from  5  =24.6  deg  for  Case  27b  to  =  26.4  deg  for  Case  29b. 

The  pore  water  pressure  acting  on  Section  C-C  remained  unchanged  since  the 
elevation  of  the  water  table  within  the  backfill  was  the  same. 

Forces  on  Section  A-A 

81.  The  changes  in  the  vertical  and  lateral  effective  forces  on  Sec¬ 
tion  A-A  were  very  nearly  the  same  as  for  Section  C-C:  Fv  increased  by 

5  percent  and  F^  decreased  by  4  percent;  the  mobilized  angle  of  friction 
increased  from  16.8  to  18.3  deg;  the  value  of  the  EPI  incre  3.ssd  f rent  0.21  for 
Case  27b  to  0.27  for  Case  29b,  indicating  a  slightly  larger  level  of  soil- 
structure  interaction. 

Forces  on  Section  H-H 

82.  As  for  Section  B-B,  the  stresses  developed  on  vertical  Section  H-H 
are  independent  of  the  lateral  movement  of  the  wall.  The  at-rest  earth 
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pressure  coefficient  (KQ)  is  equal  to  0.61  for  the  three  cases.  Fv  is  equal 
to  zero,  resulting  in  values  of  Kv  and  Sm  equal  to  zero. 

Forces  on  Section  E-E 

83.  The  variation  in  the  forces  acting  on  the  vertical  section  through 

the  toe  of  the  wall  (Section  E-E)  with  pool  elevation  are  shown  in  Figure  12. 
The  decrease  in  the  pool  elevation  is  accompanied  by  a  decrease  in  the  hydro¬ 
static  water  force  Fw  acting  on  Section  E-E.  As  Fw  decreases,  the  shear 
and  lateral  effective  forces  (F  and  Fp  acting  on  the  front  of  the  wall 
change  very  little.  This  is  due  to  the  relatively  small  increase  in  wall 
movements.  As  the  pool  of  water  was  lowered,  increased  from  0.60  to 

0.67,  Kv  decreased  from  0.13  to  0.12,  and  Sm  decreased  from  12.6  to 

10.1  deg.  However,  the  magnitude  of  the  shear  force  (T)  acting  on  the  base  of 
;he  wall  increased  considerably. 

Conclusions  on  the  effects  of 
water  in  front  of  and  behind  the  wall 

84.  The  pool  of  water  in  front  of  the  wall  provides  a  stabilizing 
influence  on  the  wall.  A  decrease  in  the  margin  of  safety  against  wall  insta¬ 
bility  occurs  when  the  pool  is  lowered  When  the  pool  is  lowered,  a  corre¬ 
sponding  contribution  to  the  resisting  force  is  lost,  which  is  then  provided 
by  increased  shear  resistance  on  the  base  of  the  wall  and,  to  a  small  extent, 
by  the  toe- fill.  The  exact  proportion  of  the  redistribution  of  the  resisting 
force  between  che  base  shear  force  and  the  toe -fill  depends  upon  the  movements 
of  the  wall.  The  larger  the  wall  movements  during  the  lowering  of  the  pod, 
the  larger  the  contribution  of  toe -fill.  In  this  series  of  analyses  most 

the  redistributed  resisting  force  was  due  to  increased  shear  force  on  the  i. 
of  the  wall,  because  the  wall  movements  were  small. 

Effect  Of  Wall  Geometry 


85.  The  effect  of  wall  geometry  was  investigated  using  the  results  of 
two  sets  of  backfill  placement  analyses;  the  first  set  of  analyses  dealing 
with  the  effect  of  wall  shape  and  the  second  set  of  analyses  dealing  with  the 
effect  of  having  a  stepped  or  a  planar  back  on  the  wall.  The  influence  of 
wall  shape  is  determined  from  a  comparison  of  the  results  of  the  backfill 
placement  analyses  of  the  three  structures  shown  in  Figure  4  (31a) ,  Figure  1 
(base  case,  17a),  and  Figure  3  (30a).  These  walls  differ  by  the  width  of 
their  crest  (T)  varying  from  a  minimum  value  corresponding  to  0.11B,  to  0.5B, 
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and  l.OB.  In  the  second  set  of  analyses,  the  results  of  the  planar  back  wall 
shown  in  Figure  4  (31a)  is  compared  to  the  results  of  the  stepped  back  wall 
shown  in  Figure  5  (32a).  No  water  table  was  simulated  in  the  analyses,  and 
the  standard  set  of  material  parameters  were  used.  The  results  are  summarized 
in  Table  4. 

Finite  element  meshes 

86.  The  finite  element  meshes  used  for  the  analyses  are  shown  in  Fig¬ 
ure  7  (base  case)  and  Figures  13-15.  The  features  of  the  three  meshes  are 
identical,  with  the  exceptions  of  the  regions  near  the  walls. 

87.  The  mesh  used  to  model  the  rectangular  wall  shown  in  Figure  13  has 
a  total  of  273  elements;  72  elements  model  the  wall,  117  model  the  backfill, 
and  40  model  the  toe-fill.  The  remaining  elements  model  the  interfaces  be¬ 
tween  the  foundation  and  the  overlying  structure  and  the  fills. 

88.  Figure  14  shows  the  245 -element  mesh  used  to  model  the  sloping  face 
wall  with  T/B  =0.11  ;  44  elements  model  the  wall,  117  model  the  backfill,  40 
the  toe -fill,  and  the  remaining  elements  model  the  interface  between  the  three 
different  material  regions. 

89.  A  total  of  274  elements  were  used  to  model  the  stepped  face  wall 
shown  in  Figure  15;  45  elements  model  the  wall,  153  model  the  backfill,  and  40 
model  the  toe -fill.  The  remaining  elements  model  all  material  interfaces.  No 
interface  elements  were  included  along  the  back  of  the  wall. 

Forces  on  Section  B-B 

90.  The  magnitude  of  the  resultant  lateral  forces  acting  on  Section  B-B 
(F^)  are  equal  for  the  three  walls  with  planar  backs  and  varying  crest  widths, 
due  to  the  fact  that  the  same  soil  properties  were  used  in  the  three  analyses. 
The  lateral  earth  pressure  coefficient  (K^)  corresponds  to  KQ  and  equals 
0.51. 

Forces  on  the  wall 

91.  Both  the  magnitudes  of  the  normal  forces  (N)  acting  along  the  base 
and  the  values  of  xR  (the  distance  from  the  toe  to  N)  increased  with  in¬ 
creasing  crest  width.  N  and  xn  vary  from  minimum  values  of  104,000  lb  and 
1.7  ft  for  T/B  =0.11  ,  to  maximum  values  of  121,800  lb  and  4.7  ft  for 

T/B  =1.0  .  The  increase  in  the  values  of  N  and  xn  with  increasing  crest 
width  reflects  the  increase  in  the  mass  of  the  wall,  accompanied  by  a  decrease 
in  the  magnitudes  of  the  shear  and  normal  forces  acting  on  the  back  of  the 
wall. 
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Summary  of  Results  for  Backfill  Placement  Analyses- -Additional  Structures 
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Section  C-C,  Soil-to-Concrete 

Interface,  Heel  Backfill  Wedge  Section  A-A,  Vertical  Section  Through  Heel 
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Soil  -  Concrete 


Figure  13.  Finite  element  mesh  used  to  model  a  rectangular  hypothetical 
structure  with  additional  backfill  beyond  the  toe  (Case  30a) 
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Soil  -Concrete 


Figure  14.  Finite  element  mesh  used  to  model  a  sloping  face  hypothetical 
structure  with  additional  backfill  beyond  the  toe  (Case  31a) 


scale 

30' 


Figure  15.  Finite  element  mesh  used  to  model  a  stepped- face  hypothetical 
structure  with  additional  backfill  beyond  the  toe  (Case  32a) 

92.  The  effective  base  contact  area  (Bg/B)  increased  from  0.36  to  0.86 
as  the  crest  width,  or  equivalently  the  mass  of  the  wall,  increased.  The 
mobilized  friction  angle  along  the  base  (5^)  decreased  from  20.7  to  15.5  deg. 
The  lateral  deformation  at  the  crest  of  the  wall  decreased  by  a  factor  of 
almost  3-1/2  (ux/H  from  0.0001  to  0.000029).  As  expected,  these  results  indi¬ 
cate  that  the  margin  of  wall  stability  increases  with  increasing  wall  mass. 
Forces  on  Section  C-C 

93.  The  magnitude  of  the  resultant  shear  force  on  Section  C-C  (Fg)  de¬ 
creased  bj  19  percent  and  the  resultant  normal  force  (Fn)  decreased  by  29  per¬ 
cent  as  the  back  of  the  wall  approaches  vertical.  The  point  of  action  of  the 
normal  force  was  constant  at  14  ft  above  the  heel.  Figure  16  shows  the  in¬ 
crease  in  mobilized  friction  angle  (§m)  with  increasing  crest  widths,  T  .  Sm 
increased  from  23.2  to  26.3  deg  (13  percent)  as  Section  C-C  approaches 
vertical. 

Forces  on  Section  A-A 

94.  The  vertical  resultant  force  on  Section  A-A  (F  )  and  the  corre¬ 
sponding  shear  stress  coefficient  (Kv)  increased  by  a  factor  of  2-1/3  as  the 
back  of  the  wall  approaches  vertical  (Figure  16) .  The  lateral  earth  pressure 
coefficient  (K^)  decreased  from  =  0.48  (Case  31a)  to  =  0.42 
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Mobilized  Angle  of 


(Case  30a) ,  corresponding  to  a  12 -percent  decrease  in  the  value  of  the 
resultant  force  .  The  mobilized  angle  of  friction  for  Section  A-A 
increased  from  10.9  to  26.3  deg. 

95.  The  distributions  of  shear  stresses  along  vertical  planes  within 
the  backfill  and  corresponding  values  of  Ky  ,  decrease  in  magnitude  with  in¬ 
creasing  distances  from  the  heel  of  the  wall,  as  shown  in  Figure  17  for  the 
rectangular  wall  (Case  30a) .  At  a  distance  of  AO  ft  the  shear  stresses  are 
nearly  equal  to  zero.  An  interesting  observation  is  that  the  height  of  the 
wall  is  also  equal  to  40  ft. 

96.  Figure  18  shows  the  variation  in  the  values  of  Ky  with  distance 

from  the  heel  of  the  wall  for  the  three  walls.  It  is  observed  that  the  value 
of  Kv  becomes  larger  as  the  back  of  the  wall  approaches  vertical.  in¬ 

creased  from  a  value  equal  to  0.09  (T/B  -  0.11),  to  -  0.14  (T/B  -  0.5), 
and  ultimately  Kv  -  0.21  (T/B  -  1.0).  In  addition,  the  distance  between  the 
heel  of  the  wall  and  the  vertical  plane  on  which  the  shear  stress  is  zero 
(equivalently  Kv  <=  0)  and  increases  as  the  width  at  the  crest  of  the  wall 
increases . 

Forces  on  Sections  H-H  and  E-E 

97.  The  lateral  earth  pressure  coefficient  (Kv)  for  Section  H-H  is  0.6 
and  equal  to  Kq  ,  since  the  base  case  material  parameters  were  used  in  the 
three  analyses. 

98.  The  shear  forces  on  Section  E-E  (F  )  are  nearly  a  constant  value 

of  3,300  lb,  with  a  corresponding  value  of  Kv  equal  to  0.16,  for  the  three 
analyses.  The  lateral  earth  pressures  showed  a  slight  variation,  with  the 
values  of  between  0.54  and  0.56,  indicating  that  the  lateral  force  (F^) 

is  more  sensitive  to  wall  displacements  than  Fv  on  Section  E-E. 

Conclusions  on  the 
effect  of  wall  shape 

99.  The  margin  of  wall  stability  increased  with  increasing  wall  mass, 
as  expected.  The  shear  force  on  the  vertical  plane  through  the  heel  of  the 
wall  increased  and  the  normal  force  decreased,  both  changes  contributing  to  an 
increase  in  wall  stability. 

Effect  of  a  stepped  back  on  the  wall 

100.  The  results  of  the  backfill  placement  analyses  for  the  wall  with  a 
stepped  back  (Case  32a)  are  summarized  in  Table  4.  These  results  may  be  com¬ 
pared  to  the  results  for  Case  31a,  in  which  the  crest  width  was  the  same 
(1.8  ft),  and  the  back  of  the  wall  was  planar. 
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Distance  From  Heel  (ft) 

Figure  18.  Variations  of  vertical  earth  pressure  coefficient 
with  distance  from  the  heel  of  the  wall  for  walls  of 

different  face  slopes 

101.  Due  to  the  fact  that  same  soil  properties  were  assigned  to  the 
backfill,  the  magnitude  of  the  resultant  lateral  forces  on  Section  B-B  are  the 
same  for  both  analyses  with  a  corresponding  lateral  earth  pressure  coefficient 
(Kq)  equal  to  0.51. 

102.  The  effective  base  contact  area  increased  by  a  factor  of  2  and  the 
maximum  compressive  stress  decreased  by  a  factor  of  2  for  the  stepped  back 
wall  (32a)  as  compared  to  the  results  of  the  planar  back  wall  (31a) .  The 
mobilized  friction  angle  along  the  base  decreased  from  20.7  to  18.5  deg,  re¬ 
flecting  both  the  decrease  in  the  value  of  base  shear  and  the  increase  in  the 
value  of  base  normal  force.  These  observations  indicate  that  the  presence  of 
an  irregular  back  on  the  wall  results  in  a  more  stable  monolith. 

103.  The  presence  of  a  stepped  back  results  in  a  decrease  in  the  mobil¬ 
ized  angle  of  friction  acting  on  Section  C-C  from  23.2  to  19.6  deg,  reflecting 
the  decrease  in  the  values  of  both  the  shear  and  normal  forces. 

104.  The  presence  of  a  stepped  back  complicates  the  stress  distribu¬ 
tions  within  the  backfill  compared  to  that  for  a  wall  with  a  planar  back. 


This  is  believed  to  be  not  only  an  inherent  feature  introduced  by  the  presence 
of  irregular  wall  geometry  adjacent  to  the  backfill  but  also  the  result  of  the 
analytical  boundary  condition  specified  along  the  interface.  The  nodal  points 
along  the  stepped  interface  were  shared  by  both  the  wall  and  the  backfill, 
introducing  a  kinematic  constraint  since  the  predominant  orientation  of  back¬ 
fill  displacemerts  are  downward  while  those  of  the  wall  are  lateral.  This  is 
in  contrast  to  the  planar  back  wall  for  which  interface  elements  were  placed 
along  all  faces  of  the  wall,  allowing  for  relative  displacements  along  the 
interface  between  the  backfill  and  the  wall.  For  the  stepped  back  wall  (32a) 
Kv  is  within  the  range  of  values  from  0.12  to  0.17  and  -  0.45  along  Sec¬ 
tion  A-A,  while  for  the  planar  back  wall  (31a)  Kv  -  0.09  and  -  0.48  .  A 
more  complete  discussion  of  the  details  regarding  the  accuracy  of  the  results 
is  given  in  Appendix  A. 

105.  The  results  of  these  analyses  show  that  a  wall  with  a  stepped  back 

has  a  larger  margin  of  stability  than  a  corresponding  wall  with  a  planar  back. 
The  presence  of  a  stepped  back  results  in  larger  values  of  Kv  and  smaller 
values  of  along  Section  A-A. 

Summary 

106.  In  the  backfill  placement  analyses  discussed  in  the  first  phase 
study  (Ebeling  et  al.  in  preparation),  a  downdrag  force  was  observed  on  the 
back  of  the  wall,  resulting  from  the  settlement  of  the  fill  under  its  own 
weight.  The  magnitude  of  the  downdrag  force  was  nearly  the  same  for  a  range 
of  values  of  Poisson's  ratio,  unit  weight  of  the  backfill,  base  interface 
shear  stiffness,  and  depth  of  toe -fill.  Fv  is  conveniently  described  in 
terms  of  a  vertical  shear  stress  coefficient  (Ky) ;  the  value  of  Kv  was  ob¬ 
served  to  be  within  the  narrow  range  from  0.13  to  0.14  for  a  number  of  cases 
in  which  the  values  of  these  parameters  were  varied. 

10/ .  In  the  backfill  placement  analyses  discussed  in  paragraphs  13 
through  69,  variations  in  the  magnitudes  of  the  stiffness  of  the  backfill  or 
interface  shear  stiffness  resulted  in  larger  variations  in  the  value  of  Kv  , 
ranging  from  0.09  to  0.15.  In  the  evaluation  of  walls  v?ith  different  geome¬ 
tries,  the  magnitude  of  Kv  varied  with  the  inclination  of  the  face  of  the 
wall,  increasing  in  value  from  0.09  to  0.21  as  the  back  of  the  wall  approached 
vertical.  In  addition,  a  stepped  face  wall  had  a  larger  value  of  Kv  than  a 
planar  face  wall. 

108.  An  increase  in  the  stiffness  of  the  backfill,  or  a  reduction  in 
the  values  of  the  soil-to-concrete  interface  shear  stiffnesses  resulted  in  a 
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slight  increase  in  the  lateral  earth  pressures  acting  on  the  wall.  As  the 
back  of  the  wall  approached  vertical,  the  lateral  earth  pressures  decreased. 

As  expected,  the  margin  of  wall  stability  increased  with  increasing  wall  mass. 

109.  A  pool  of  water  in  front  of  the  wall  provides  a  stabilizing  influ¬ 
ence.  A  decrease  in  the  pool  elevation  results  in  reduction  of  the  resisting 
force  and  an  increase  in  base  shear  and  the  force  exerted  by  the  toe-fill. 

The  distribution  of  the  resisting  force  between  base  shear  and  toe-fill  resis¬ 
tance  depends  on  the  magnitude  of  the  movements  that  result  from  lowering  the 
water  level  in  front  of  the  wall. 
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PART  III:  EVALUATION  OF  THE  STABILITY  OF  RETAINING 
STRUCTURES  USING  CONVENTIONAL  EQUILIBRIUM  METHODS 

110.  Part  III  describes  evaluations  of  the  stability  of  earth- retaining 
structures  founded  on  rock  using  conventional  equilibrium  methods.  In  Part  VI 
of  the  Ebeling  et  al.  (in  preparation)  initial  phase  of  study,  a  comparison 
between  the  results  of  CEA  and  FEA  was  made  for  the  base  case  structure.  The 
CEA  and  FEA  results  were  extended,  as  described  in  these  paragraphs,  to 
include  consideration  of  the  behavior  of  the  backfill  and  toe-fill.  The 
interdependency  between  the  mobilized  shear  resistances  on  the  base  of  the 
wall  and  the  earth  pressures  within  the  soil  fills  were  investigated.  The 
results  were  summarized  using  conventional  concepts  of  the  factor  of  safety 
against  shear  failure  along  the  base  of  the  wall  (sliding)  and  within  the 
backfill  (referred  to  as  shear  factors),  related  to  the  magnitudes  of  the 
forces  acting  on  the  wall. 

111.  An  additional  set  of  analyses  are  described  in  which  the  defor¬ 
mations  of  a  retaining  structure  during  backfilling  are  sufficient  to  break 
the  bond  between  the  base  of  the  wall  and  its  rock  foundation.  The  influence 
of  the  base  shear  properties  on  the  base  of  the  wall  was  examined. 

Factor  of  Safety  and  Shear  Factors 


112.  Earth- retaining  structures  of  the  type  discussed  in  this  report 
can  be  considered  in  three  parts:  the  wall,  the  backfill,  and  the  toe-fill. 
The  backfill  applies  the  driving  force  while  the  wall  and  the  toe- fill  provide 
the  resisting  forces.  In  a  CEA,  assumptions  are  made  regarding  the  magnitude 
of  the  forces  applied  by  the  backfill,  because  the  conditions  of  equilibrium 
are  insufficient  for  determining  the  magnitude  of  the  load  without  assump¬ 
tions.  The  magnitude  of  these  assumed  earth  loads  correspond  to  some  level  of 
mobilized  shear  resistance  within  the  backfill.  The  magnitude  of  the  earth 
loads  can  be  described  using  a  shear  factor  (SF).  For  each  fill  region  there 
exists  an  SF,  defined  as  the  ratio  of  the  ultimate  shear  force  along  a  poten¬ 
tial  slip  plane  divided  by  the  shear  force  required  for  equilibrium. 

113.  When  earth  loads  on  a  retaining  wall  are  assumed  and  the  factor  of 


safety  against  sliding  (F£>base 


)  is  computed,  corresponding  values  of  the  SF's 


for  the  backfill  and  toe -fill  can  be  inferred  from  the  magnitude  of  the  earth 
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loads.  This  interdependence  between  these  factors  has  been  examined  using  CEA 
procedures . 

114.  The  base  case  structure  without  toe-fill  is  shown  in  Figure  19. 


8  =  16  ft. 

Figure  19.  Base  case  hypothetical  structure 
without  toe-fill 

The  structure  is  16  ft  wide  at  the  base,  8  ft  wide  at  the  crest,  and  retains 
40-ft  of  backfill.  In  the  Case  13a  CEA,  summarized  in  Table  7,  Part  VI  of  the 
Ebeling  et  al.  (in  preparation)  first  phase  study,  the  at- rest  earth  pressure 
force,  corresponding  to  K  -  0.51  ,  was  applied  to  Section  A-A,  the  vertical 
plane  through  the  heel  of  the  wall.  Using  the  CEA  forces  on  the  base  of  the 
structure,  the  mobilized  base  friction  angle  was  30.5  deg,  which  is  equal  to 
the  assumed  angle  of  sliding  resistance  for  the  interface.  That  is,  the  fac¬ 
tor  of  safety  against  sliding  along  the  base  (FSbase)  was  equal  to  1.0. 

115.  The  assumptions  regarding  the  magnitude  of  the  forces  acting  on 
Section  A-A  reflect  a  corresponding  value  of  the  mobilized  shear  resistance 

within  the  backfill,  characterized  by  the  SF  (SF,  ,  For  the  cohesion- 

J  backfill 

less  backfill  of  Case  13a,  the  SF  is  given  by 


‘backfill  tan  <f>' 
m 
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where 


4>'  -  effective  angle  of  internal  friction 
-  mobilized  angle  of  internal  friction 

116.  A  limiting  case  occurs  when  the  strength  of  the  soil  is  fully 

mobilized  (^F^  “  1-0)  i  which  coincides  with  the  development  of  active 

earth  pressures  within  the  backfill. 

117.  The  mobilized  angle  of  friction  (</>^)  for  this  case,  corresponding 

to  -  0.51  and  Kv  -  0  ,  is  equal  to  IS  deg.  This  value  of  is 

determined  using  the  Rankine  relationship 


(8) 


This  relationship  is  only  valid  when  Ky  is  equal  to  zero.  In  cases  where 
Kv  is  not  zero,  another  equilibrium  relationship  would  be  used.  Coulomb's 
relationships  for  and  Kv  ,  discussed  by  Kezdi  (1975)  ,  can  be  used  when 

Ky  has  a  value  greater  than  zero. 

118.  With  <f> '  equal  to  39  deg  for  the  backfill,  the  SFbackfill  » 
defined  by  Equation  7,  is  equal  to  2.36.  It  can  be  observed  that  although  the 
wall  is  on  the  verge  of  sliding  with  the  assumed  earth  load,  the  backfill  has 
a  significant  margin  of  safety  against  shear  failure. 

119.  In  the  previously  described  CEA,  the  assumed  set  of  forces  acting 

on  Section  A-A  (K^  "0.51  and  Ky  =  0)  corresponded  to  ^backfill  “2.36 

and  FS,  "1.0  ,  Variations  in  the  value  assumed  for  Ku  ,  between  K 
base  h  o 

(0.51)  and  (0.23),  results  in  the  relationship  between  the  values  of 

^backfill  aTK*  ^base  shown  in  Figure  20.  It  can  be  observed  that  when  the 
strength  of  the  backfill  is  fully  mobilized  ^Fbackfill  ”  ancl  ^h  “  ^a  ’ 

the  wall  has  an  ample  margin  of  safety  against  sliding  (FS^ase  =  2.24).  The 

values  of  ^Fbackfill  an<*  ^base  wou^  be  ecJual  (1-47)  if  was  equal 

to  0.35. 


120.  When  assumptions  are  made  regarding  the  magnitude  of  the  forces 
acting  on  a  wall,  a  level  of  mobilized  shear  resistance  within  the  backfill  is 
implied,  and  this  value  may  be  characterized  by  a  shear  factor.  This  shear 
factor  may  be  determined  using  conventional  earth  pressure  theories.  There  is 
thus  an  interdependence  between  the  computed  factor  of  safety  against  failure 
in  shear  along  the  base  of  the  wall  and  the  shear  factor  for  the  backfill. 
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sz 

X. 


Figure  20.  Variations  of  the  lateral  earth 
pressure  coefficient  and  the  factor  of 
safety  against  base  shear  failure  with 
shear  factor  of  the  backfill 
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Factor  of  Safety.  Shear  Factors.  Load  Factor,  and 
Resistance  Ratio 


121.  In  FEA's,  the  loading  exerted  by  the  backfill  on  the  wall  is 
determined  by  the  interaction  between  the  wall  and  the  backfill  during  place¬ 
ment  of  the  backfill  behind  the  wall.  The  finite  element  method  of  analysis 
differs  from  CEA  in  that  it  is  not  necessary  to  estimate  or  assume  the  magni¬ 
tude  of  the  loads  imposed  on  the  structure  by  the  backfill. 

122.  The  results  from  three  backfill  placement  analyses  of  the  base 
case  structure  were  summarized  in  Table  7  of  the  first  phase  study  (Ebeling  et 
al.  in  preparation).  No  water  pressures  were  represented  in  the  analyses. 

The  effect  of  soil  fill  at  the  toe  of  the  wall  was  investigated  using  the 
results  of  Case  13a  (without  toe-fill)  and  Case  17a  (with  17.8  ft  of  toe- 
fill).  In  Case  19a  the  shear  stiffness  assigned  to  the  interface  along  the 
base  of  the  wall  was  reduced  from  its  standard  value,  resulting  in  increased 
wall  movements,  reduced  backfill  loads,  and  greater  toe-fill  resistance 
forces . 

123.  The  effective  force  acting  on  Section  A-A  (called  F  )  is  the 
resultant  of  the  forces  Fv  and  .  The  resultant  effective  earth  pres- 
sure  coefficient  (K  )  can  be  defined  as 


* 


K 


* 

F 


rH  - 

Jo  % dh 


(9) 


Using  the  relationships  presented  in  Part  II, 


K 


can  be  expressed  as 


* 


K 


(10) 


Table  5  shows  that  the  values  of  K  for  the  backfill  are  nearly  equal  to 
,  due  to  the  small  contribution  of  the  shear  force  acting  on  Section  A-A. 
The  mobilized  angles  of  internal  friction  ( <f> ^)  within  the  backfill  were 
computed  using  Figure  6  of  the  Department  of  the  Navy  NAVFAC  DM7.2  Manual 
(1982).  Coulomb's  relationships,  discussed  by  Kezdi  (1975),  may  also  be  used, 
and  they  result  in  the  same  values  of  <f> ^  .  Figure  21  shows  the  resulting 
variation  of  the  shear  factor  with  the  mobilized  angle  of  friction  of  the 


57 


Table  5 

Shear  Factor.  Load  Factor,  and  Resistance  Ratio 
for  Backfill  and  Toe-Fill 

a.  Shear  and  Load  Factors  for  Backfill 


Case 

Section  A-A 

Backfill 

K 

V 

*h 

K* 

8 

m 

deg 

rm 

deg 

Shear 

Factor 

^backfill 

K 

a 

Load 

Factor 

13a 

0.138 

0.453 

0.474 

16.9 

19 

2.35 

0.23 

2.06 

17a 

0.137 

0.454 

0.474 

16.8 

19 

2.35 

0.23 

2.06 

19a 

0.138 

0.337 

0.364 

22.2 

25 

1.74 

0.23 

1.58 

Shear  Factor  and  Resistance  Ratio  for  Toe-Fill 


Section 

E-E 

Toe -Fill 

Case 

K 

V 

“h 

K* 

5 

m 

deg 

<y 

m 

deg 

Shear 

Factor 

SF 

toe 

K 

P 

Resistance 

Ratio 

13a 

-- 

-- 

-- 

-- 

-- 

-- 

-- 

17a 

0.156 

0.548 

0.569 

15.9 

14 

3.25 

14.24 

0.04 

19a 

0.069 

1.446 

1.448 

2.7 

12 

3.81 

14.24 

0.10 

Note:  Shear  factor  =  tan  <f>'/ tan  4> ^  where  <f>'  ■=  39  deg 
Load  factor  =  K*/K 

a 

Resistance  ratio  ■=■  K*/K 

P 

--  indicates  not  applicable 


backfill.  It  is  observed  that  as  the  mobilized  angle  of  friction  increases, 
the  resultant  force  applied  to  the  wall  decreases;  that  is,  the  wall  does  less 
of  the  work  resisting  the  applied  force  as  more  of  the  work  is  done  by  the 
backfill. 

124.  For  Cases  13a  and  17a,  <f>^  was  equal  to  19  deg,  and  for 

Case  19a,  ^  was  equal  to  25  deg  as  shown  in  Table  5a.  Using  Equation  7, 
with  the  effective  angle  of  internal  friction  equal  to  39  deg,  the  resulting 
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Figure  21.  Variations  of  shear  factor  and 
load  factor  with  mobilized  angle  of 
friction,  backfill 


shear  factors  for  the  backfills  (SF^  were  found  to  be  equal  to  2.35 

and  1.74  as  shown  in  this  table.  From  Table  7  (Ebeling  et  al.  in  prepara¬ 
tion),  FS^asg  “  1-29  for  Case  13a,  FS^ase  -  1.75  for  Case  17a,  and 
^bace  ”  for  Case  19a.  For  the  base  case  analyses  (13a  and  17a),  the 

margin  of  safety  is  greater  for  the  backfill  than  the  base,  where  the  reverse 
is  observed  for  Case  19a. 

125.  Another  index  that  can  be  used  to  characterize  the  magnitude  of 
the  applied  load  is  referred  to  as  the  load  factor,  defined  as  the  ratio  of 
the  resultant  earth  pressure  coefficient  (K  )  divided  by  the  active  earth 
pressure  coefficient  (K  =  0.23).  The  variation  in  the  load  factor  with 

a 

mobilized  angle  of  friction  for  the  backfill  is  shown  in  Figure  21.  The 
larger  the  value  of  <f>^  ,  the  smaller  the  load  applied  to  the  wall  (K  ) .  The 
limiting  value  of  load  factor  equal  to  unity  corresponds  to  the  full  mobiliza¬ 
tion  of  the  shear  resistance  within  the  backfill.  In  this  condition,  the 
earth  pressures  are  equal  to  their  active  values.  For  Cases  13a  and  17a,  the 
load  factor  was  equal  to  2.06,  and  for  Case  19a,  the  load  factor  was  equal 
to  1.58. 


126.  The  presence  of  17.8  ft  of  toe- fill  (Cases  17a  and  19a)  provides 
additional  resistance  to  the  loads  applied  by  the  backfill,  and  decreases  the 
contribution  required  from  the  wall  to  the  total  resistance.  The  proportion 
of  total  resisting  force  provided  by  the  toe -fill  depends  on  the  magnitude  of 
the  wall  deformations;  increased  wall  deformations  result  in  larger  contribu¬ 
tions  from  the  toe-fill. 

127.  Characterizing  the  magnitude  of  the  force  applied  by  toe- fill  is 

more  involved  than  was  the  case  for  the  backfill.  The  value  of  <j>'  is 

dependent  upon  the  orientation  of  the  shear  force  developed  along  the  slip 

plane  (T*)  as  shown  in  Figure  22.  When  the  shear  force  (T*)  counteracts  the 

effect  of  the  gravity  force,  the  relationship  for  the  shear  factor  (sFtoe)  to 

the  left  of  <b'  •=*  0  is  valid. 
rm 

128.  When  T  acts  in  the  direction  consistent  with  W  ,  the  relation- 

ship  to  the  right  of  <f>^  =  0  is  valid.  In  this  case,  the  value  of  K  is 
greater  than  1.0.  The  equilibrium  relationship  used  is  that  associated  with 
the  development  of  passive  pressures,  which  is  sensitive  to  the  value  of  the 
mobilized  angle  of  friction  on  Section  E-E  (Sm). 

129.  The  reversal  in  the  direction  of  T  occurs  when  K  -  1.0  ,  for 
which  T  -  0  and  <f>^  =  0  ,  resulting  in  an  infinite  value  for  the  shear 
factor.  Thus,  the  shear  factor  has  two  limitations  as  an  index;  it  is 
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Resistance 
Ratio,  K/Kp  1.0 


0.2 


discontinuous ,  and  each  value  represents  two  different  levels  of  wall  loads , 
corresponding  to  two  different  directions  of  the  shear  stress  on  a  potential 
slip  plane. 

130.  Another  index  that  can  be  used  for  the  toe- fill,  which  overcomes 
the  limitations  of  the  shear  factor,  is  referred  to  as  the  resistance  ratio. 
The  resistance  ratio  is  defined  as  the  ratio  of  the  resultant  effective  earth 
pressure  coefficient  (K  )  for  Section  E-E  divided  by  the  passive  earth  pres¬ 
sure  coefficient  (Kp  -  14.24  for  S/<j>'  -  0.8).  Figure  22  shows  the  variation 
in  the  resistance  ratio  with  mobilized  angle  of  friction  for  the  toe -fill.  It 
can  be  seen  that,  unlike  the  shear  factor,  the  resistance  ratio  is  continuous 
and  .single  valued.  The  resistance  ratio  has  a  minimum  value  of 

Ka/K  -  0.02  ,  when  the  earth  pressures  are  equal  to  their  active  value 

•k  * 

(K  -  Ka) ,  and  a  maximum  value  of  unity,  when  the  earth  pressures  are  equal  to 
their  passive  value  (K  -  Kp). 

131.  As  noted  previously,  the  active  and  passive  states  of  stress  both 
correspond  to  full  mobilization  of  the  shear  resistance  within  the  backfill 
but  with  T*  acting  in  opposite  directions.  In  Case  17a,  the  resistance 
ratio  of  the  toe-fill  is  equal  to  0.04.  The  resistance  ratio  increases  to  a 
value  of  0.10  as  a  result  of  the  increased  wall  deformations  in  Case  19a 
(Table  5).  Although  there  was  a  two  and  one-half  fold  increase  in  the  resis¬ 
tance  ratio  in  Case  19a,  90  percent  of  the  resisting  force  provided  by  toe- 
fill  still  has  not  been  mobilized. 

132.  In  the  base  case  analysis  (17a),  the  driving  force  of  the  backfill 


was  characterized  by  SFbackfiH  “  ^*35  ’  t^ie  Foai*  factor  equal  to  2.06, 

while  the  resisting  forces  of  the  base  and  the  toe -fill  were  characterized  by 


FSbase  “  1.75  and  the  resistance  ratio  of  the  toe-fill  was  equal  to  0.04. 

With  the  increased  wall  deformations  for  Case  19a,  SF,  ,  =1.74 

backfill 

(26  percent  decrease),  the  load  factor  was  1.58  (23  percent  decrease), 


FS,  =11.6  ,  and  the  resistance  ratio  was  0.10. 
base 


It  may  be  observed  that 


increased  wall  deformations  resulted  in  an  increase  in  the  degree  of  soil 


structure  interaction,  exemplified  by: 


a. 


b. 


Mobilization  of  a  greater  portion  of  the  available  shear 
resistance  within  the  backfill,  as  evinced  by  the  decrease  in 


the  value  of 


^backfill 


A  reduction  in  the  magnitude  of  the  load  applied  by  the  back¬ 
fill  to  the  wall,  characterized  by  a  decrease  in  the  value  of 
the  load  factor. 
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c.  An  increase  in  the  resistance  provided  by  the  toe- fill,  with  a 
corresponding  increase  in  the  value  of  the  resistance  ratio. 

d.  The  required  contribution  of  shear  on  the  wall  to  the  total 

resisting  force  was  reduced,  as  exemplified  by  the  increase  in 

the  value  of  the  FS, 

base 

Loss  of  Bond  Between  Wall  and  Foundation 


133.  This  section  describes  a  set  of  three  equilibrium  analyses  con¬ 
sidering  a  hypothetical  case  in  which  the  deformations  of  a  retaining  struc¬ 
ture  during  backfilling  are  sufficient  to  break  the  bond  along  the  interface 
between  the  base  of  the  wall  and  the  rock  foundation.  The  interdependence 
between  the  forces  acting  on  the  wall,  the  ultimate  shear  resistance  along  the 
base  of  the  wall,  and  the  resulting  deformations  were  examined  in  the 
analyses . 

134.  The  structure  shown  in  Figure  23  retains  backfill  20  ft  in  height, 


T 


Very  Stiff  Rock 


E, 


Figure  23.  Hypothetical  structure  used 
in  example  of  loss  of  bond  along  the 
base  of  the  wall 
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is  3.5  ft  thick,  12  ft  wide  at  its  base,  and  is  buttressed  by  3.5  ft  of  toe- 
fill.  The  forces  acting  on  the  structure  are: 

W  -  weight  of  the  retaining  structure 

-  14,963  lb 

E,  -  earth  load  resulting  from  the  placement  of  backfill 

E2  -  resisting  force  due  to  the  presence  of  toe- fill 
N  -  normal  force  along  the  base  of  the  wall 

-  W 

-  shear  load  along  the  base  of  the  wall 

135.  The  vertical  shear  forces  acting  on  the  back  and  the  front  of  the 
wall  were  assumed  to  be  equal  to  zero,  which  from  previous  equilibrium  analy¬ 
ses  (i.e.  Part  VI  (Ebeling  et  al.  in  preparation))  is  known  to  be  a  conserva¬ 
tive  assumption.  No  water  pressures  were  considered  in  the  analyses. 

Material  parameters 

136.  The  standard  set  of  material  parameters  were  assigned  to  the  clean 
granular  backfill  (Dr  -  75  percent) ;  the  unit  weight  of  the  soil  is  equal  to 
135  pcf  and  the  effective  angle  of  internal  friction  is  39  deg.  The  strength 
of  the  concrete -to -rock  interface  was  modeled  using  Mohr-Coulomb  strength 
parameters.  In  a  limited  number  of  direct  shear  tests  conducted  by  the 

US  Army  Corps  of  Engineers  on  the  interfaces  between  concrete  and  competent 
rock,  described  by  Benson  et  al.  (1987),  deformations  on  the  order  of  0.01  in. 
resulted  in  the  loss  of  bond  between  the  materials.  Prior  to  loss  of  bond 
between  the  concrete  and  rock,  the  peak  shear  resistance  of  the  interface  be¬ 
tween  the  base  of  the  wall  and  the  rock  foundation  (r  )  was  given  by 

t  =  C  +  cr  .  tan  6,  (11) 

p  a  n  bp 


where 

Ca  =  bond,  or  adhesion,  between  the  concrete  and  the  rock 
=  2,000  psf 

CTn  “  normal  stress  along  the  interface 
^bp  ”  Pea^  angle  of  internal  friction 

137.  The  three  analyses  differed  by  the  values  of  the  angle  of  internal 
friction  assigned  to  the  interface;  5^  =  30  deg  for  Case  1,  5^  “  15  deg 

for  Case  2,  and  6^  "  10  deg  for  Case  3. 
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138.  For  lateral  wall  deformations  greater  than  0.01  in.,  the  loss  of 
bond  along  the  base  of  the  wall  was  characterized  by  a  value  of  C  equal  to 

a 

zero  and  a  residual  shear  strength  (r  )  given  as 


n 


•  tan  8 


br 


(12) 


where  6^r  is  the  residual  angle  of  internal  friction.  The  values  of 
were  equal  to  the  values  of  5^  for  all  three  cases. 

and  E2  for 

complete  backfill  and  toe -fill 

139.  The  magnitude  of  the  lateral  forces  E^  and  E2  ,  resulting  from 

the  placement  of  fill  behind  and  in  front  of  the  wall,  may  be  expressed  using 

the  lateral  earth  pressure  coefficient  .  The  variation  of  the  lateral 

earth  pressure  coefficient  with  wall  movement  is  shown  in  Figure  24.  The 

average  yield  of  the  soil  is  defined  as  the  average  lateral  deformation  of  the 

backfill  along  the  wall-to-soii  interface  ((u  )  )  divided  by  the  height  of 

x  ave 

the  backfill  (H*) .  This  term  was  first  introduced  by  Terzaghi  (1936)  to 
explain  the  consistency  between  the  measured  lateral  load  versus  deformation 
behavior  of  the  sand  backfills  in  the  1934  tests  (Terzaghi  1934)  using  trans¬ 
lating  and  rotating  walls.  It  is  important  to  recogrize  that  the  relationship 
in  Figure  24  models  the  variation  of  with  average  yield  based  on  the 

final  fill  height  (H*) . 

140.  When  the  average  displacement  along  the  soil-wall  interface 

((uv)  )  is  equal  to  zero,  the  values  of  the  at-rest  earth  pressure  coeffi- 

x  ave 

cients  (Kq)  are  equal  to  0.5  for  the  backfill  and  0.6  for  the  toe- fill.  As 

the  wall  displaces  the  toe-fill,  idealized  at  the  upper  left  in  Figure  24, 

increases  from  a  value  of  Kq  =  0 . 6  ,  as  shown  in  this  chart  at  the  bottom  of 

this  figure.  When  (u  )  was  equal  to  4.2  in.  [  (u  )  /H*  =■  0.1],  the  earth 

x  ave  x  ave 

pressures  attained  their  passive  values  (K^  *=>4.4)  in  the  toe -fill. 

141.  As  the  wall  moved  away  from  the  backfill,  shown  at  the  upper  right 

in  Figure  24,  K-n  decreased  from  a  value  of  Ko  =  0.5  ,  as  shown  by  this 

chart  at  the  bottom  of  the  figure.  When  (ux)gyp  was  equal  to  0.96  in. 

((u  )  /H*  *=  0.004),  the  earth  pressures  were  equal  to  their  active  values 

x  ave 

(Ka  -  0.23)  within  the  backfill.  It  can  be  seen  that  the  values  of  average 
yield  required  to  attain  active  pressures  within  the  backfill  are  an  order  of 
magnitude  smaller  than  those  required  to  attain  passive  pressures. 
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Figure  24.  Variation  of  the  lateral  earth 
pressure  coefficient  with  average  yield 


142.  When  the  variation  in  Kh  of  Figure  24  is  normalized  by  its 
extreme  values,  the  development  of  the  limiting  states  of  pressures  (Kp  and 
Ka)  may  be  compared  with  those  measured  in  tests  involving  retaining  walls 
(ierzaghi  1934,  Johnson  1953)  or  computed  in  previous  analytical  studies 
(Clough  and  Duncan  1971),  as  shown  in  Figures  25  and  26. 

143.  In  Figure  25  the  reduction  in  with  average  yield  away  from 
the  backfill  is  shown.  For  the  backfill,  a  normalized  value  of  zero 
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Variation  of  normalized  lateral  earth  pressure  coe 
average  yield  of  the  wall  away  from  the  backfill 


Figure  26.  Variation  of  normalized  lateral  earth  pressure  coefficient  with 
average  yield  of  the  wall  displacing  fill 


corresponds  to  Kq  (0.50)  and  a  normalized  value  of  unity  corresponds  to 
(Kh)min  "  (0.23).  It  can  be  seen  that  the  shape  of  the  normalized  curve 

in  Figure  25  is  nonlinear.  A  significant  portion  of  the  reduction  in 
occurred  during  the  first  stages  of  wall  movement  away  from  the  backfill.  The 
shape  of  the  curve  used  in  the  analyses  was  consistent  with  the  three  tests 
conducted  by  Terzaghi  (1934)  using  a  wall  retaining  approximately  5  ft  of 
loose  or  dense  sand,  and  with  the  analytical  studies  of  Clough  and  Duncan 
(1971) ,  who  used  the  finite  element  method  to  analyze  a  wall  retaining  10  ft 
of  loose  sand. 

144.  In  Figure  26,  the  increase  in  with  movement  of  the  wall 

toward  the  toe -fill  is  shown.  For  the  toe -fill,  a  normalized  displacement 

value  of  zero  corresponds  to  Kq  (0.6)  and  a  normalized  displacement  value  of 

unity  corresponds  to  (K,  )  **  K  **  4.4  .  The  shape  of  the  normalized  curve 

n  max  P 

for  the  toe-fill  was  more  nearly  linear  than  was  the  case  for  the  backfill. 

It  can  be  observed  that  the  relationship  of  Figure  24  is  consistent  with  the 
shape  of  the  curves  for  the  six  experimental  tests  conducted  by  Johnson  (1953) 
using  a  wall  retaining  2  ft  of  dense  sand. 

Values  of  and  during  placement  of 
soil  fill  in  the  finite  element  analyses 

145.  This  section  describes  the  development  of  a  relationship  between 
the  height  of  the  fill,  the  applied  forces,  and  the  displacements  for  the 
wall.  It  is  based  on  the  results  of  the  backfill  placement  analyses  using  the 
finite  element  method  and  therefore  models  the  behavior  of  a  wall  when  there 
is  no  loss  of  bond  and  no  sliding  of  the  wall  along  the  concrete- to -rock 
interface. 

146.  Several  backfill  placement  analyses  of  retaining  walls  were  dis¬ 
cussed  in  Part  VI  of  the  first  phase  study  (Ebeling  et  al.  in  preparation)  and 
in  Part  II  of  this  report.  The  variations  of  earth  pressure  forces  with  wall 
movement  are  shown  in  Figure  27  for  the  40-ft-high  wall  of  Case  13a  (no  toe- 
fill)  and  Case  17a  (17.8  ft  of  toe-fill).  The  nine  points  represent  the  re¬ 
sults  after  placement  of  each  of  the  4.44-ft-thick  soil  lifts.  It  can  be  seen 

that  the  relationship  between  F,  and  (uv)„„  (Figure  27)  is  nonlinear. 

n  x  ave 

147.  After  completion  of  backfilling  for  Case  13a  (no  toe-fill), 

Fh  ■=  48,922  lb  (K^  =  0.45)  and  (ux)ave  -  14  x  10  ^  ft.  The  eight  intermedi¬ 
ate  values  of  F,  and  (u  )  were  normalized  by  these  final  values,  as 
n  x  oV6 

shown  in  Figure  28. 
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Figure  27.  Variation  of  lateral  force  acting 
on  plane  A-A  with  wall  displacement 

148.  For  the  base  case  structure  with  17.8  ft  of  toe-fill  (Case  17a), 

the  final  values  of  F,  -  49,008  lb  (K,  =■  0.45)  and  (u  )  =>  11  x  10"4  ft. 

The  normalized  intermediate  values  of  Fh  ,  (uy)avp  ,  and  height  of  fill  dur¬ 
ing  the  nine  stages  of  backfilling  are  shown  in  Figure  29. 

149.  Case  17a  differs  from  Case  13a  by  the  placement  of  17.8  ft  of  toe- 
fill.  During  the  first  four  stages  of  backfilling,  4.44-ft-thick  lifts  of 
fill  were  placed  in  front  of  as  well  as  behind  the  wall,  resulting  in  no 
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Figure  28.  Variation  of  normalized  lateral  force  acting 
on  plane  A-A  with  normalized  height  of  backfill  and 
wall  displacements,  without  toe-fill 
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Figure  29.  Variation  of  normalized  lateral  force  acting 
on  plane  A-A  with  normalized  height  of  backfill  and 
wall  displacements,  toe-fill 


lateral  displacements  until  the  height  of  backfill  exceeded  17.8  ft.  Due  to 
the  presence  of  the  toe -fill,  the  final  average  lateral  displacements  were 
smaller  by  22  percent. 

150.  The  normalized  results  for  the  intermediate  stages  of  backfilling 
for  Cases  19a  and  30a  are  shown  in  Figure  29.  In  Case  19a,  a  reduced  shear 
stiffness  was  assigned  to  the  concrete-to-rock  interface  for  the  base  case 
wall.  Compared  to  Case  17a,  the  backfill  load  on  Section  A-A  was  reduced  and 
the  wall  movements  were  increased. 

151.  The  rectangular  wall  of  Case  30a  retains  backfill  40  ft  in  height 
and  is  buttressed  by  17.8  ft  of  toe- fill.  For  this  wall,  both  the  lateral 
force  on  Section  A-A  and  the  average  displacement  were  smaller  than  those  of 
the  base  case.  It  can  be  seen  in  Figure  29  that  there  is  good  agreement  be¬ 
tween  the  normalized  results  for  the  three  cases,  which  is  quite  interesting 
since  the  differences  between  the  magnitudes  of  the  final  lateral  forces  and 
concrete-to-soil  interface  displacements  were  so  large. 

152.  The  normalized  values  of  the  lateral  forces  acting  on  Section  A-A 
and  average  wall  movements ,  are  shown  in  Figure  30  for  the  four  analyses .  The 
lateral  forces  are  characterized  by  the  ratio  K^/Kq  .  The  magnitude  of  the 
lateral  force  acting  on  Section  A-A  decreases  with  increasing  wall  movement 
and  is  lower  when  the  back  of  the  wall  is  vertical. 

153.  The  results  for  Cases  17a,  19a,  and  30a  v/ere  used  to  estimate  the 
increase  in  lateral  displacements  for  walls  without  toe -fill  by  extrapolation 
of  the  computed  results  for  walls  with  toe-fill,  as  shown  in  Figure  27  for 
Case  17a.  The  procedure  resulted  in  a  predicted  increase  in  the  total  dis¬ 
placements  of  7  x  10  5  ft.  It  can  be  seen  in  Figure  28  that  the  resulting 
normalized  relationship  for  Case  17a  is  in  good  agreement  with  that  of 

Case  13a,  the  base  case  wall  without  toe-fill.  This  is  due  in  part  to  the 
smaller  increments  of  wall  movements  occurring  during  the  first  stages  of 
backfill  placement,  as  compared  to  the  latter  stages.  The  extrapolation 
procedure  was  also  applied  to  the  results  from  Cases  19a  and  30a,  as  shown  in 
Figure  28.  The  variation  in  the  results  for  the  normalized  forces  (during 
backfilling)  versus  normalized  displacements  was  larger  than  that  of  Fig¬ 
ure  29,  but  the  variation  in  normalized  forces  with  height  of  backfill  was 
about  the  same. 

154.  Using  the  normalized  relationships  of  Figures  28  and  29,  the  rela¬ 
tionships  for  the  wall  shown  in  Figure  22  (htQe/H  =  0.175)  was  developed,  as 
shown  in  Figure  31.  To  use  Figure  31,  the  values  of  F^  on  Section  A-A  and 
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Figure  30.  Normalized  lateral  earth  pressure  coefficient 
versus  average  yield  of  wall  40  ft  in  height 

(ux)ave  >  after  completion  of  backfilling,  are  required.  These  were  deter¬ 
mined  using  Figure  30.  For  a  40-ft-high  wall  with  a  vertical  back  and  an 
intermediate  concrete- to-rock  interface  shear  stiffness,  the  average  yield 
(after  backfilling)  was  equal  to  40  x  10 ’ 5  and  K^/KQ  is  equal  to  0.66.  In 
Part  V  of  the  first  phase  of  study  (Ebeling  et  al.  in  preparation),  it  was 
shown  that  the  displacement  of  a  wall  varied  with  the  square  of  its  height. 
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The  average  yield  would  therefore  be  approximately  20  x  10  for  a  wall  20  ft 
in  height,  with  the  final  (ux)ave  equal  to  0.048  in. 

155.  The  precise  value  of  the  factor  between  displacements  for  like 

structures  of  different  heights  is  not  known  and  thus  needs  to  be  evaluated  in 
future  studies.  Since  KQ  was  equal  to  0.5  in  the  analyses,  was  equal 

to  8,950  lb  after  completion  of  20  ft  of  backfilling.  Intermediate  values  of 

and  (ux)ave  during  backfilling  are  determined  using  the  normalized  rela¬ 
tionships  shown  in  Figure  31. 

Loss  of  bond  along  base  of 
wall- -equilibrium  analyses 

156.  Three  backfill  placement  analyses  were  performed  for  the  retaining 
structure  shown  in  Figure  23,  considering  possible  deformations  during 
backfilling  exceeding  0.01  in.,  which  would  be  sufficient  to  break  the  bond 
between  the  base  of  the  wall  and  the  rock  foundation.  The  interdependence 
between  the  forces  a«_t?ng  on  the  front  and  back  of  the  wall  (E^  and  £2),  the 
ultimate  shear  resistance  along  the  base  of  the  wall  (Tb)  and  the  resulting 
deformations  can  be  demonstrated  by  these  equilibrium  analyses.  The  three 
analyses  differed  by  the  value  of  the  angle  of  friction  assigned  to  the  inter¬ 
face  (6bp  and  Sbr) ;  the  values  of  5bp  and  6br  were  equal  to  30  deg  for 
Case  1,  15  deg  for  Case  2,  and  10  deg  for  Case  3. 

Case  1  equilibrium  analysis 

157.  The  forces  acting  on  the  wall  (Figure  23)  at  the  end  of  backfill¬ 
ing  and  the  subsequent  lateral  wall  displacements  are  shown  in  Figure  32  for 
the  Case  1  analysis  (5bp  *»  5br  -  30  deg).  E^  is  the  earth  load  resulting 
from  the  placement  of  backfill,  E2  is  the  passive  resistance  from  the  toe- 
fill,  and  Tb  is  the  shear  load  acting  along  the  base  of  the  wall. 

158.  The  dashed  relationship  labeled  (Ei)compiete  backfill  describes 
the  variation  in  E^  ,  after  completion  of  backfilling  with  the  average 
lateral  wall  displacements,  as  given  in  Figure  24.  The  earth  load  (E^)  must 
intersect  this  relationship  at  the  completion  of  backfilling,  which  is  labeled 
Stage  5  in  the  figures. 

159.  The  ultimate  resistance  of  the  wall,  labeled  "shear  and  passive 
resistance"  in  Figure  32,  is  equal  to  the  sum  of  the  ultimate  shear  force 
developed  along  the  base  of  the  wall  and  E2  .  This  relationship  is  discon¬ 
tinuous  at  a  lateral  displacement  equal  to  0.01  in.,  the  displacement  that 
corresponds  to  the  loss  in  bond  between  the  wall  and  the  foundation.  For 
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Figure  32.  Variation  of  forces  acting  on  the  wall  with  deformation,  Case 


displacements  less  than  0.01  in.,  the  ultimate  shear  resistance  is  given  by 
Equation  11  with  C&  -  2,000  psf  (1  tsf)  and  5^  -  30  deg.  At  0.01  in.  of 
displacement,  the  adnesion  was  lost  (Ca  -  0),  and  the  ultimate  shear  resis¬ 
tance  is  given  by  Equation  12  with  Sbr  -  30  deg. 

160.  For  a  bonded  base,  the  value  of  E^  for  any  level  of  backfill  is 
determined  using  Figure  31,  with  E^  ,  after  completion  of  backfilling,  equal 
to  8,950  lb  and  the  final  (ux)ave  equal  to  0.048  in.  (corresponding  to  a 
final  average  yield  equal  to  20  x  10  "*).  This  relationship  is  valid  until  the 
loss  of  bond  occurs  ((ux)ave  “0.01  in.)  and  under  certain  conditions,  to  be 
discussed  subsequently,  is  valid  after  the  displacements  exceed  0.01  in. 

161.  The  backfill  placement  equilibrium  analyses  involved  as  many  as 
five  stages.  The  first  stage  corresponded  to  the  simultaneous  placement  of 
backfill  and  toe-fill  to  a  height  equal  to  ^toe  (3.5  ft),  with  no  lateral 
displacement.  The  value  of  E^  was  equal  to  537  lb,  using  Figure  31,  and  E2 
was  equal  to  496  lb,  using  Figure  24.  By  equilibrium  of  forces  in  the  hori¬ 
zontal  direction,  the  shear  load  along  the  base  of  the  wall  is  equal  to  the 
difference  between  E^  and  Eg  ,  41  lb.  The  ultimate  shear  and  passive 
resistance  (33,135  lb)  was  much  larger  than  the  required  resistance  (T^  +  Eg 

-  537  lb),  and  as  expected,  no  shear  failure  occurred  along  the  base  of  the 
wall  at  this  stage  of  backfilling. 

162.  As  backfilling  behind  the  wall  proceeded,  the  wall  moved  away  from 
the  backfill,  displacing  the  toe -fill.  The  value  of  E^  increased  with  in¬ 
creasing  height  of  backfill  (Figure  31) ,  and  the  value  of  Eg  (Figure  24) 
increased  with  increasing  wall  displacements.  With  the  increase  in  the  magni¬ 
tude  of  E^  greater  than  that  of  Eg  ,  an  increase  in  the  shear  load  on  the 
base  (T^)  resulted.  When  the  height  of  backfill  was  equal  to  14.8  ft 

(h/H  -  0.74),  the  displacement  of  the  wall  was  equal  to  0.01  in.  (Figure  31), 
resulting  in  the  loss  of  bond  along  its  base.  At  this  stage  of  backfilling, 
labeled  Stage  2  in  the  Figure  32,  E-^  -  5,102  lb,  Eg  -  580  lb,  and 

Tb  =  4,522  lb. 

163.  The  transition  in  the  ultimate  shear  resistance  from  the  peak 
value  to  the  residual  value  is  shown  in  Figure  33  for  an  average  normal  stress 
equal  to  1,250  psf  along  the  base  of  the  wall.  The  ultimate  shear  stress- 
deformation  behavior  of  the  interface  was  assumed  to  be  brittle,  as  observed 
in  the  direct  shear  test  results  reported  by  Benson  et  al.  (1987).  The  ulti¬ 
mate  shear  resistance  was  reduced  from  a  value  of  2,722  to  722  psf,  due  to  the 
loss  of  the  bond  (C&  «*  2,000  psf).  This  resulted  in  a  decrease  in  the 
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ultimate  shear  and  passive  resistance  from  33,219  to  9,219  lb,  as  shown  in 
Figure  32.  Since  the  ultimate  shear  and  passive  resistance  (9,219  lb)  is 
larger  than  the  required  resistance  (T^  +  E^  -  5,102  lb),  shear  failure  did 
not  occur  on  the  base  of  the  wall. 

164.  After  completion  of  backfilling  (Stage  5),  the  final  deformation 
of  the  wall  was  equal  to  0,048  in.,  with  E-^  -  8,950  lb,  E2  “  900  lb,  and 
-  8,050  lb.  The  ultimate  shear  and  passive  resistance  (9,539  lb)  was 


Figure  33.  Available  shear  resistance  during  loss  of 
bond  between  the  wall  and  the  foundation 
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larger  than  the  required  resistance  (T^  +  E2  -  8,950  lb),  thus  no  shear  fail¬ 
ure  occurred  along  the  base  of  the  wall  at  any  stage  of  backfilling.  Since  no 
shear  failure  occurred  during  the  simulated  placement  of  backfill,  the  rela¬ 
tionships  of  Figure  31  remained  valid  throughout  the  analysis 
Case  2  equilibrium  analysis 

165.  In  the  second  bac^-Llll  placement  analysis  of  the  wall  shown  in 

Figure  23,  all  the  material  parameters  were  the  same  as  those  used  in  Case  1, 
with  the  exception  that  the  values  of  the  peak  and  residual  base  friction 
angles  were  equal  to  15  deg.  The  relationships  between  the  forces  acting  on 

the  wall  (E^  ,  E2  ,  and  T^)  and  the  average  lateral  displacement  of  the  wall 

are  shown  in  Figure  34.  The  Figure  24  relationship  between 

^complete  backfill  and  ^Vave  >  described  in  tbe  Case  1  analysis,  was 
also  used  in  developing  the  relationships  snown  in  Figure  34.  The  ultimate 

shear  and  passive  resistance  of  the  wall  was  computed  using  the  same  procedure 
as  for  Case  1,  but  with  5,  and  Sbr  equal  to  15  deg, 

166.  The  forces  acting  on  the  wall  and  the  displacements  of  the  wall 
after  Stage  1,  the  completion  of  filling  in  front  of  and  behind  the  wall  to 
h  -  h  e  (3.5  ft),  were  the  same  as  in  the  previous  analysis;  =  537  lb, 

E2  -  496  lb,  T^  -  41  lb,  and  (ux)ave  "  0  The  ultimate  shear  and  passive 
resistance  (28,505  lb)  was  much  large c  than  the  required  resistance  (T^  +  E2 
=  537  lb) ,  thus  no  shear  failure  occurred  along  the  base  of  the  wall  at  this 
stage  of  backfilling. 

167.  At  Stage  2,  when  the  height  of  backfill  was  equal  to  14.8  ft  (h/H 

-  0.7*/,  E^  =  5,102  lb,  E2  -  580  lb,  Tb  =  4,522  lb,  and  (t^,)  =0.01  in. 

They  are  the  same  set  of  values  as  computed  in  the  Case  1  analysis  since  the 
same  relationships  for  E^  (Figure  31)  ai.d  for  Eg  (Figure  24)  were  used  in 
both  analyses.  At  this  value  of  lateral  displacement,  the  base  of  the  wall 
loses  its  bond  with  the  foundation,  resulting  in  a  transition  in  the  ultimate 
shear  resistance  from  peak  (C&  =  2,000  psf  and  5bp  =  15  deg)  to  residual 
values  (Ca  =  0  and  6br  =  15  deg) .  A  corresponding  decrease  in  the  ultimate 
shear  and  passive  resistance  from  28,589  to  4,589  lb  occurred,  as  shown  in 
Figure  34  (labeled  N  •  tan  Sbr,  +  Eg  in  the  figure  with  the  expanded  scale). 
Since  the  value  for  the  ultimate  shear  and  passive  resistance  (4,589  lb)  was 
less  than  the  required  resistance  (5,102  lv) ,  a  shear  failure  along  the  base 
of  the  wall  resulted.  Subsequently,  the  wall  ioved  away  from  the  backfill, 
displacing  the  toe-fill.  The  additional  displ "cements  resulted  in  a  decrease 
in  the  value  of  E^  and  an  increase  in  K.e  /al*<.  of  Eg  .  The  movement 
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34.  Variation  of  forces  acting  on  the  wall  with  deformation,  Case 


terminated  when  equilibrium  was  attained  between  the  ultimate  shear  and  pas¬ 
sive  resistance  and  the  required  resistance  (T^  +  .  Equilibrium  was  rees¬ 

tablished  after  a  45 -percent  increase  in  the  lateral  displacement  of  the  wall, 
to  0.0145  in.  At  this  stage,  labeled  Stage  4  in  Figure  34,  the  value  of 
was  4,602  lb.  Thus,  equilibrium  among  the  forces  acting  on  the  wall  was 
attained  by  reducing  the  earth  load  and  increasing  the  passive  resistance 
provided  by  the  toe-fill,  through  the  lateral  translation  of  the  wall. 

168.  Since  the  loss  of  bond  along  the  base  resulted  in  a  shear  failure 

and  subsequent  lateral  wall  movements  to  reestablish  equilibrium  among  the 
forces,  the  relationships  between  the  height  of  backfill,  the  lateral  force, 
and  the  height  of  backfill  and  lateral  displacements  shown  in  Figure  31  were 
no  longer  valid.  The  increase  in  the  value  of  E2  with  lateral  displacements 
controlled  the  behavior  of  the  wall  when  backfilling  resumed,  since  the  shear 
load  (T^)  was  equal  to  the  ultimate  shear  resistance  (N  •  tan  5^  -  4,009  lb) 

along  the  base  of  the  wall.  During  this  phase  of  backfilling,  the  earth  load 
(E^)  was  equal  to  the  shear  and  passive  resistance.  The  final  lateral  dis¬ 
placement  of  „he  wall  was  determined  by  the  intersection  of  the  shear  and 
passive  resistance  relationship  with  the  relationship  for 

(El>complete  backfill'  The  value  of  <ux>final  was  ec* ual  to  °'55  in' 

(Stage  5),  “  6,409  lb  and  E2  *=»  2,400  lb. 

169.  In  this  case,  where  a  reduced  angle  of  internal  friction 

(5  =  15  deg)  was  assigned  to  the  base,  a  loss  in  bond  along  the  concrete- to- 
rock  interface  resulted  in  shear  failure  during  backfilling,  with  a  lateral 

translation  of  the  wall  required  to  reestablish  equilibrium  among  the  loading 

and  resisting  forces.  The  behavior  of  the  wall  at  subsequent  stages  of  back¬ 
filling  was  controlled  by  the  residual  shear  strength  along  the  base  and  pas¬ 
sive  force  provided  by  the  toe -fill.  The  final  displacement  of  the  wall  in 
this  case  was  11  times  as  large  as  that  for  Case  1. 

Case  3  equilibrium  analysis 

170.  The  results  of  the  third  backfill  placement  analysis  of  the  wall 
(Figure  23)  are  shown  in  Figure  35.  This  equilibrium  analysis  differed  from 
the  previous  Case  1  and  Case  2  analyses  by  the  value  of  8^  and  5^r  ,  each 
equal  to  10  deg.  From  Figure  24,  the  relationship  for  (E1)complete  backfill 
was  obtained,  and  the  ultimate  shear  and  passive  resistance  was  computed  as 
described  in  the  previous  two  analyses . 

171.  At  Stage  1,  after  backfilling  to  3.5  ft  (h  =  htoe) ,  E^  =  537  lb, 

E0  -  496  lb,  Tu  -  41  lb,  and  (u  )  -  0  .  The  ultimate  shear  and  passive 
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35.  Variation  of  forces  acting  on  the  wall  with  deformation.  Case 


resistance  was  equal  to  27,134  lb,  while  the  required  resistance  (Tb  +  E2)  was 
equal  to  537  lb. 

172.  After  14.8  ft  of  backfilling  was  completed  and  the  lateral  dis¬ 
placement  of  the  wall  was  equal  to  0.01  in.  (Stage  2),  -  5,102  lb, 

E2  ■»  580  lb,  Tb  -  4,522  lb,  as  in  Cases  1  and  2.  The  ultimate  shear  and  pas¬ 
sive  resistance  decreased  from  27,218  to  3,218  lb,  with  the  loss  of  the  bond 
along  the  base  of  the  wall.  Since  the  value  of  the  required  resistance,  equal 
in  value  to  E^  ,  was  larger  than  the  sum  of  N  •  tan  $br  and  E2  ,  as  shown 
in  the  upper  figure  with  the  expanded  scale,  a  shear  failure  along  the  base  of 
the  wall  resulted.  As  in  Case  2,  the  wall  moved  away  from  the  backfill,  dis¬ 
placing  the  toe-fill;  thereby  decreasing  the  value  of  E-^  and  increasing  the 
value  of  E2  . 

173.  When  (ux)ave  “  0-028  in.,  the  earth  pressures  were  equal  to  their 
active  values  (K  =  0.23),  and  En  attained  a  minimum  value  of  3,401  lb 
(Stage  3).  The  wall  continued  to  move  laterally  because  the  required  resis¬ 
tance,  equal  to  3,401  lb  (same  value  as  E^) ,  was  less  than  the  ultimate  shear 
and  passive  resistance,  equal  to  3,370  lb.  Since  equilibrium  was  not 
achieved,  the  wall  continued  to  move,  resulting  in  a  further  increase  in  the 
value  of  E2  ,  while  the  values  of  E^  (3,401  lb)  and  Tb  (-  N*tan  $br 

-  2,638  lb)  were  constant.  Equilibrium  was  reestablished  at  (ux)ave 

-  0.032  in.  (Stage  4),  when  the  required  resistance  (■*  E^)  was  the  same  magni¬ 
tude  as  the  ultimate  shear  and  passive  resistance,  a  value  of  3,401  lb.  Thus 
equilibrium  was  attained  through  the  translation  of  the  wall  in  two  stages: 
first,  a  reduction  in  the  earth  load  to  active  pressures,  accompanied  by  an 
increase  in  the  value  of  the  resistance  provided  by  the  toe-fill,  and 
secondly,  a  continued  increase  in  the  value  of  the  toe-fill  resistance,  with  a 
constant  earth  load.  During  the  wall  movements,  the  shear  load  acting  along 
the  base  of  the  wall  was  equal  to  the  ultimate  residual  shear  resistance. 

174.  The  relationships  between  height  of  backfill,  wall  load,  and  wall 
displacement  of  Figure  31  were  no  longer  valid  since  a  shear  failure  occurred 
along  the  base  of  the  wall,  as  in  the  Case  2  analysis.  The  increase  in  the 
value  of  the  toe-fill  resistance  (E2)  with  displacements  controlled  the  subse¬ 
quent  behavior  of  the  wall  when  backfilling  was  resumed.  During  the  final 
stages  of  backfilling,  the  shear  load  along  the  base  of  the  wall  (Tb)  was 
equal  to  the  ultimate  shear  resistance  (2,638  lb),  and  E^  was  equal  to  the 
sum  of  Tb  and  E2  ,  as  shown  in  Figure  35.  The  final  lateral  wall 
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displacement  was  equal  to  3.3  in.  (Stage  5).  It  is  observed  that  at  (ux)ave 

-  0.96  in.  f (u  )  /H  -  0.004],  the  earth  load  attained  its  minimum  value  of 
6,210  lb,  which  corresponds  to  the  active  earth  pressures  (K  -  0.23).  During 

a 

placement  ot  the  last  5.2  ft  of  backfill,  from  Stage  4  to  Stage  5,  the  back¬ 
fill  was  in  a  state  of  active  failure,  and  the  base  shear  load  was  equal  to 
the  ultimate  residual  shear  resistance. 

175.  A  reduced  ultimate  shear  resistance  along  the  base  (5^  and 

-  10  deg)  resulted  in  larger  lateral  wall  translations  than  the  previous 

analyses.  In  fact,  the  displacement  was  69  times  larger  than  that  of  Case  1 
(S^p  and  5^r  -  30  deg) .  The  larger  lateral  displacements  for  Case  3  were 
required  to  reestablish  equilibrium  among  the  forces  acting  on  the  wall  during 
backfilling  when  a  loss  of  bond  and  subsequent  shear  failure  occurred  along 
the  base.  Like  Case  2,  the  subsequent  behavior  of  the  wall  was  controlled  by 
the  residual  shear  strength  along  the  base  and  the  passive  force  provided  by 
the  toe-fill.  The  final  displacement  of  the  wall  was  larger  by  a  factor  of 

about  6.0  than  that  of  Case  2  (S^p  anc*  ^br  "  ^  >  because  the  displace¬ 

ments  during  backfilling  were  of  a  sufficient  magnitude  so  as  to  achieve  a 
minimum  value  for  the  earth  load,  corresponding  to  active  earth  pressures  (Kfl 

-  0.23),  when  the  height  of  backfill  exceeded  14.8  ft. 
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PART  IV:  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Summary 

176.  The  results  discussed  in  this  second  phase  describe  the  continued 
evaluations  of  the  accuracy  and  reliability  of  CEA  of  earth- retaining  struc¬ 
tures  founded  on  rock. 

177.  As  a  result  of  several  simplifying  assumptions,  CEA  of  gravity 
earth- retaining  structures  neglect  the  true  process  of  soil-structure  interac¬ 
tion  and  tend  to  produce  highly  conservative  results.  It  is  believed  that 
improvements  based  on  soil -structure  interaction  studies  like  those  described 
in  this  report  are  likely  to  lead  to  more  realistic  and  less  conservative 
predictions  of  response.  This  assertion  is  supported  by  the  fact  that  some 
existing  gravity  walls  under  Corps  of  Engineers'  jurisdiction,  e.g.  navigation 
lock  walls,  do  not  satisfy  current  stability  criteria  based  on  CEA,  while  at 
the  same  time  they  do  not  exhibit  signs  of  distress  or  substandard 
performance . 

178.  The  studies  performed  under  this  research  contract  are  divided 
into  first  and  second  phases.  The  first  phase  report  (Ebeling  et  al.  in  prep¬ 
aration)  defined  the  analytical  requirements  for  this  category  of  structure 
and  resulted  in  the  development  of  new  finite  element  procedures,  while  this 
report  explored  a  number  of  variables  on  the  stability  of  the  wall.  The 
results  of  the  initial  phase  were  discussed  by  Ebeling  et  al.  (in  preparation) 
while  the  results  of  this  study,  the  second  phase,  are  included  in  both 
reports . 

179.  The  gravity  earth -retaining  structures  evaluated  in  this  study 
were  consistent  with  the  characteristics  of  Corps  of  Engineers  retaining 
structures  as  described  in  the  initial  report.  In  addition,  an  "L"  shaped 
earth-retaining  structure  was  used  in  the  studies  of  the  loss  of  bond  and 
possible  slippage  along  the  base  of  the  wall  during  backfilling. 

180.  The  results  from  FEA  and  equilibrium  analyses  utilizing  the 
results  from  select  FEA  were  used  in  the  evaluations.  These  studies  were 
performed  using  "backfill  placement"  analysis  procedures  in  which  the  place¬ 
ment  of  backfill  behind  the  walls  was  simulated,  layer  by  layer.  The  finite 
element  program  used  in  this  series  of  analyses  has  the  capability  to  model 
the  development  of  a  gap  between  the  base  of  the  wall  and  its  rock  foundation, 
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and  this  phenomenon  was  analyzed  in  Case  4  as  described  in  the  initial  phase 
study. 

181.  In  the  parametric  studies  discussed  in  this  volume,  it  was  ob¬ 
served  that  increases  in  stiffness  of  the  backfill  or  a  reduction  in  the  val¬ 
ues  of  soil-to-concrete  interface  shear  stiffnesses  resulted  in  a  slight 
increase  in  the  lateral  earth  pressures  acting  on  the  wall.  Variations  in  the 
magnitudes  of  the  stiffness  of  the  backfill  or  the  interface  shear  stiffnesses 
resulted  in  larger  variations  in  the  magnitude  of  the  shear  force  acting  on  a 
vertical  plane  passing  through  the  heel  of  the  wall  (Section  A-A) .  In  the 
initial  analyses,  it  was  observed  that  vertical  shear  stress  coefficient  (Ky) 
was  within  the  narrow  range  of  0.13  to  0.14  for  variations  in  the  magnitudes 
of  Poisson's  ratio,  the  unit  weight  of  the  backfill,  and  base  interface  shear 
stiffness.  Studies  discussed  in  this  report  indicated  that: 

a.  The  value  of  increased  from  0.12  to  0.15  as  the  stiffness 

of  the  backfill  was  reduced  over  the  range  of  values  of  stiff¬ 
ness  representative  of  clean  granular  backfills. 

b.  Ky  decreased  in  value  from  0.14  to  0.09  as  the  values  of  the 
concrete-to-soil  and  rock-to-soil  interface  shear  stiffnesses 
were  reduced. 

182.  A  reduction  in  the  magnitude  of  the  stabilizing  shear  force  acting 
on  the  wall  results  in  a  reduction  of  the  actual  margin  of  wall  safety. 

183.  As  the  slope  of  the  back  of  the  wall  increased,  the  magnitude  of 
K  increased  from  a  value  equal  to  0.09  to  a  value  equal  to  0.21  for  a  wall 
with  a  vertical  back  face.  As  expected,  it  was  found  that  the  margin  of  wall 
stability  increased  with  increasing  wall  mass.  In  an  additional  set  of  analy¬ 
ses,  it  was  observed  that  a  stepped  face  wall  has  a  greater  margin  of  stabil¬ 
ity  than  a  smooth  face  wall  due  to  a  larger  value  of  Kv  and  slightly  smaller 
earth  pressure  loads. 

184.  In  the  evaluation  of  the  base  case  structure  with  water  in  front 
of  and  behind  the  wall,  it  was  observed  that  the  water  in  front  of  the  wall 
provides  a  stabilizing  influence  ,  decrease  in  the  pool  elevation  results  in 
the  loss  of  a  corresponding  contribution  to  the  resisting  force,  which  is  then 
provided  by  shear  on  the  base  of  the  wall  and  the  toe -fill.  The  exact  propor¬ 
tion  of  the  redistribution  of  the  resisting  force  between  the  wall  and  the 

toe -fill  depends  upon  the  movements  of  the  wall.  The  greater  the  wall  move¬ 
ments,  the  larger  the  contribution  of  the  toe- fill  to  the  redistributed  re¬ 
sisting  force. 
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185.  Results  from  the  CEA  and  the  FEA  were  used  to  show  the  interdepen¬ 
dency  between  the  factors  of  safety  for  the  structure  and  the  backfill.  The 
factor  of  safety  against  a  shear  failure  within  the  backfill  was  referred  to 
as  a  shear  factor.  It  was  observed  that  the  factor  of  safety  against  sliding 
for  the  structure  and  the  shear  factor  of  the  backfill  are  not  necessarily 
equal  in  magnitude. 

186.  The  interdependence  between  wall  deformations  and  the  loading  and 
resisting  forces  was  examined  in  a  series  of  equilibrium  analyses  in  which  the 
deformation  of  the  wall  during  backfilling  was  hypothesized  to  be  sufficient 
to  break  the  bond  between  the  base  of  the  wall  and  the  rock  foundation.  It 
was  demonstrated  that  a  backfill  placement  analysis  can  be  conducted  using 
equilibrium  methods,  and  that  the  magnitude  of  the  ultimate  base  shear  resis¬ 
tance  and  the  resistance  provided  by  the  toe -fill  controls  the  displacement  of 
the  wall.  When  a  loss  of  bond  along  the  base  of  the  wall  occurs  during  back¬ 
filling,  the  wall  must  slide  along  its  base  to  reestablish  equilibrium. 

Two  actions  occur  during  sliding:  the  earth  load  is  reduced  as  a  result  of 
the  wall  movement  away  from  the  backfill,  while  the  passive  force  acting  on 
the  toe  of  the  wall  is  increased  with  the  wall  movement  toward  the  toe -fill. 
During  sliding,  the  strength  of  the  interface  along  the  base  is  at  its 
residual  value.  In  general,  the  lower  the  available  resistance  along  the 
interface  between  the  wall  and  its  foundation,  the  greater  the  resulting  wall 
displacements . 

187.  Finally,  this  study  has  shown  that  the  displacements  of  the  wall 
have  a  significant  influence  on  the  distribution  of  both  stabilizing  and 
destabilizing  forces  exerted  on  the  wall  by  the  fill,  and  on  the  base  of  the 
wall.  In  general,  as  the  wall  moves  away  from  the  backfill,  the  lateral 
forces  behind  the  wall  decrease  and  the  lateral  forces  in  front  of  the  wall 
increase.  Calculations  indicate  that  the  magnitude  of  the  wall  displacements 
required  to  develop  significant  changes  in  the  magnitude  of  the  shear  forces 
acting  on  the  base  of  the  wall  are  less  than  those  required  to  develop  active 
pressures  behind  the  wall  and  much  less  than  those  required  to  develop  passive 
pressures  in  front  of  the  wall. 


Conclusions 


188.  The  results  from  these  studies  indicate  that  the  traditional  as¬ 
sumptions  used  in  the  CEA  of  earth-retaining  structures  founded  on  rock  are 
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conservative.  One  of  the  assumptions  in  the  CEA  is  that  the  shear  force 
acting  on  the  plane  through  the  heel  of  the  wall  is  equal  to  zero.  In  the 
finite  element  studies  for  a  select  set  of  hypothetical  walls,  with  character 
istics  similar  to  those  of  Corps  of  Engineers  structures,  it  was  found  that 
the  shear  force  acting  on  a  vertical  plane  within  the  backfill,  passing 
through  the  heel  of  the  wall,  was  not  zero.  These  studies  showed  that  the 
values  of  Kv  for  the  walls  studied  varied  from  0.09  to  0.21,  depending  on 
the  geometrical  and  material  parameters  of  the  wall  and  backfill. 

Recommendations 


189.  Because  of  the  potential  importance  of  these  findings  with  regard 
to  evaluations  of  the  stability  of  gravity  retaining  structures,  it  is  recom¬ 
mended  that  the  studies  be  continued  to  investigate  two  important  questions 
which  remain  unanswered: 

a.  Will  similar  analytical  results  be  found  for  walls  founded  on 
soil? 

b.  Can  the  existence  of  shear  forces  on  the  backs  of  gravity 
walls  be  confirmed  experimentally? 
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APPENDIX  A:  INTERPRETATION  OF  THE  FINITE  ELEMENT  RESULTS 


1.  In  the  traditional  finite  element  method  of  analysis  the  primary  un¬ 
knowns  are  the  nodal  point  displacements,  and  the  secondary  unknowns  are  the 
stresses  within  the  elements.  For  each  of  the  two-dimensional  elements  used 
by  the  program  SOILSTRUCT,  stresses  are  computed  only  at  the  centers  of  the 
elements,  although  the  stress  varies  across  the  element.  This  variation  may 
be  nonlinear,  depending  upon  the  distortion  of  the  element.  However,  the 
variation  in  stress  across  an  interface  element  is  linear,  due  to  the  mathe¬ 
matical  formulation  of  this  element. 

2.  The  results  of  the  finite  element  analyses  discussed  in  this  and  the 
initial  study  (Ebeling  et  al.  in  preparation)*  are  summarized  using  the  resul¬ 
tant  forces  acting  along  tl.e  interfaces  between  the  wall  and  the  backfill,  the 
wall  and  the  rock  foundation,  and  along  vertical  planes  within  the  backfill. 
For  the  planes  passing  through  the  two-dimensional  elements,  such  as  vertical 
Section  A- A  within  the  backfill,  the  resultant  forces  were  computed  using  the 
stresses  calculated  at  the  center  of  each  of  the  soil  elements,  as  shown  in 
Figure  Al  (Method  1) .  These  computations  were  made  by  first  graphically  con¬ 
structing  shear  and  lateral  stress  distributions  along  Section  A-A,  connecting 
the  stresses  between  element  centers  as  shown  in  the  figure.  The  set  of  re¬ 
sultant  forces  (Fv  and  F^)  and  point  of  action  (hp^)  were  then  computed  from 
the  plotted  stress  distributions. 

3.  Along  the  base  of  the  wall  and  vertical  Section  E-E,  the  resulting 
stress  distributions  are  known  since  the  values  of  stresses  are  computed  at 
the  nodal  points  and  their  variation  across  the  element  is  linear.  Thus,  the 
computed  magnitudes  and  point  of  action  of  the  resulting  interface  forces  are 
exact. 

4.  An  alternative  procedure  for  computing  the  magnitude  of  the  resul¬ 
tant  forces  and  their  point  of  action  on  Section  A-A  (Method  2)  is  shown  in 
Figure  Al.  In  this  procedure  the  unknown  values  of  Fv  ,  F^  ,  and  hp^  are 
determined  using  the  three  equations  of  equilibrium  (vertical,  horizontal,  and 
moment) ,  since  the  remaining  forces  of  the  weight  and  resultant  interface 
forces  along  the  base  and  Section  E-E  are  known  both  in  magnitude  and  point  of 
action. 


*  References  cited  in  this  Appendix  are  included  in  the  References  following 
the  main  text. 
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re  Al.  Two  methods  for  computing  the  magnitude  and  point  of  action  of 
the  resultant  forces  along  plane  A-A 


5.  The  resultant  forces  acting  on  Section  A-A  were  computed  using  both 
methods  for  all  of  the  backfill  placement  analyses  discussed  herein  and  in  the 
initial  study  (Ebeling  et  al.  in  preparation),  and  the  results  are  summarized 
in  Table  Al.  It  is  observed  that  the  values  of  lateral  resultant  forces  (F^) 
and  their  points  of  action  along  Section  A-A  (h^)  are  in  agreement  for  the 
two  methods  in  all  cases.  The  values  of  the  resultant  shear  forces  (Fy)  com¬ 
puted  by  the  two  methods  differ  by  less  than  8  percent  for  all  except 

Cases  24a  and  32a.  This  corresponds  to  a  difference  in  values  of  vertical 
earth  pressure  coefficients  (Ky)  of  less  than  0.01. 

6.  Case  24a  is  a  parametric  evaluation  with  stiffness  parameters  that 
do  not  reflect  those  characteristic  of  typical  backfills.  In  Case  32a,  the 
difference  between  the  values  of  Fv  reflect  an  involved  stress  field,  due  to 
both  the  presence  of  irregular  geometry  adjacent  no  the  backfill  (i.e.  stepped 
interface)  and  an  analytical  boundary  condition  specified  along  this  inter¬ 
face.  For  this  analysis,  the  nodal  points  along  the  stepped  interface  were 
shared  by  both  the  wJl  and  the  backfill.  This  introduced  a  kinematic  con¬ 
straint  since  the  predominant  direction  of  backfill  displacements  are  down¬ 
ward,  while  those  of  the  wall  are  lateral  and  directed  away  from  the  backfill. 

7.  It  was  concluded  that  the  accuracy  of  the  results  are  dependent  not 
only  upon  the  interpretation  of  the  resulting  stress  field  but  also  the  size 
of  the  mesh.  For  Case  32a,  it  is  concluded  that  the  true  value  of  vertical 
earth  pressure  coerficient  is  within  the  range  from  0.12  to  0.17. 
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Summary  of  Results  for  Backfill  Placement  Analyses 


Cu 

w 


HONr4H00HC'Jf'^<j  f*4 

C'JCMCNCMCvJr-<VOC'4r-lr-<r-4i-«r-ICMC'JC'4C^COr-Jf-< 

o  ooooooooooooooooooo 


A4 


for  Cases  27b,  28b,  and  29b. 


APPENDIX  B:  SOILSTRUCT  (ISOTROPIC)  USER'S  MANUAL  WITH  BASE 
SEPARATION  MODEL  USING  THE  ALPHA  METHOD 


1.  IDENTIFICATION  CARD  FORMAT  (20A4) 


Column 

Variable 

Explanation 

1-80 

HED 

Analysis  identification. 

2 .  DATA 

CONTROL  CARD 

FORMAT  (1315,  D5.1,  D10.1) 

(a)  MESH  PARAMETERS 

All  nodes  and  two-dimensional  and  interface  elements  to  be  used  in  the 
analysis  must  initially  be  included  in  the  mesh;  additions  or  deletions 
are  not  allowed,  but  the  material  parameters  can  be  changed  to  make  them 
inactive.  One -dimensional  bar  elements  usually  are  added  in  subsequent 
construction  steps,  but  may  initially  be  included  in  the  mesh. 

Column 

Variable 

Explanation 

1-5 

NUMNP 

Number  of  nodal  points  (500  maximum) . 

6-10 

NUMEL 

Number  of  elements  initially  in  the  mesh, 
excluding  bar  and  beam  elements ,  but  including 
interface  elements.  NUMEL,  then,  includes  NUMJT 
(450  maximum) .  Interface  elements  should  be 
numbered  first. 

11-15 

NUMJT 

Number  of  interface  elements  (50  maximum) . 

16-20 

NUMBAR 

Number  of  bar  elements  initially  in  the  mesh 
(15,  including  elements  added  in  subsequent  con¬ 
struction  steps) . 

(b)  ANALYSIS  PARAMETERS 

Column 

Variable 

Explanation 

21-25 

NC 

Number  of  loading  and  construction  steps 
(60  maximum) . 

26-30 

NMOD 

Modulus  specification  code 

=  0  if  modulus  calculation  codes  input 
with  Loading  Information  Card. 

«=  1  if  modulus  calculation  codes  input 
with  Modulus  Calculation  Card. 
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Column 

Variable 

Exolanation 

31-35 

INIT 

Initial  stress  input  code 

“0  if  external  input  from  cards  or 
tape  included  in  input  sequence. 

=1  if  internally  generated  from 
gravity  turn-on  analysis. 

-  2  if  the  initial  stresses  are  to  be 
set  equal  to  zero  internally. 

(If  INIT  is  not  equal  to  0,  1,  or  2,  initial  stresses 
are  generated  assuming  a  horizontal  ground  surface  with 
the  water  table  at  the  ground  surface,  both  soil  and 
interface  elements  are  included  in  the  mesh.  In  this 
case,  initial  stresses  are  generated  by  the  direct 
gravity  (a  -  Kq  a  )  fcr  a  horizontal  ground  and  water 
table.  The  unit  weight  of  the  soil  is  that  of  the 
first  soil  for  which  properties  are  given.) 

36-40 

KI 

Equal  to  zero  (interface  is  automatically 
activated) .  KI  is  used  when  INIT  is  not  equal 
to  0.  If  INIT  =  0,  then  KI  can  be  left  blank. 

41-45 

IHORIZ 

Ground  surface  inclination  code  - 

■=  0  horizontal  ground  surface.  Verti¬ 
cal  stresses  computed  from  gravity 
turn-on.  Horizontal  stresses  are 
computed  assuming  K?  =  i//(l  -  v) . 

Unless  Kq  is  specified. 

=  1  sloping  ground  surface.  Vertical 
and  horizontal  stresses  are  calcu¬ 
lated  from  gravity  turn-on  analy¬ 
sis  linear  elastic  response  soil, 
i.e. ,  Kq  =  u/(l  -  i/) . 

(c)  OUTPUT 

PARAMETERS 

Column 

Variable 

Explanation 

46-50  ITRD  Analysis  printout  code 

■=0  if  initial  stresses  are  to  be 
printed. 

=2  if  initial  stresses  are  not  to  be 
printed. 

51-55  ILIST  Element  and  nodal  point  card  data 

printout  code 

=0  if  not  printed. 

=1  if  printed. 

56-60  I PUNCH  Idle. 
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Column 


Variable 


Explanation 


61-65 


ITAPE 


Initial  disc  storage  code 
-  0  if  no  disc  storage. 

”1  if  storage  of  analysis  continua¬ 
tion  data  and  displacements  from 
final  iterations  is  required. 


BASIC  PARAMETERS 


The  continuation  data  and  displacements  are 
stored  in  separate  files  that  are  set  up  on  the 
IBM  job  control  cards  (tape  9). 


The  unit  weight  of  water  and  the  atmospheric  pressure  are  included  as 
basic  parameters.  Either  English  or  SI  units  can  be  used,  but  must,  of 
course,  be  compatible  with  input  coordinate,  pressure,  and  material 
property  parameters. 


Column 


Variable 


Explanation 


66-70 


/ 1-80 


Unit  weight  of  water. 
Atmospheric  pressure. 


3.  MATERIAL  ALLOCATION  CARD  FORMAT  (615) 

All  two-dimensional  material  types  are  assigned  numbers  first,  followed  by  the 
interface  material  types.  (Note  that  bar  elements  are  not  assigned  a  material 
type  number,  but  are  identified  solely  by  their  element  number.)  If  a  number 
for  NATYP,  for  example,  is  not  required,  assign  a  value  of  zero  or  leave 


blank. 

Column 

Variable 

Explanation 

1-5 

NUMMAT 

Total  number  of  material  types,  including 
two-dimensional  soil  or  construction  material 
types  and  interface  material  types  (40  maximum) . 

6-10 

NUMSOL 

Total  number  of  material  types  excluding 

the  interface  material  types.  Thus  the  differ¬ 
ent  between  NUMMAT  and  NUMSOL  must  equal  the 
number  of  interface  material  types. 


11-15  NATYP  Material  type  number  assigned  to  two-dimensional 

elements  having  the  properties  of  air.  Usually, 
elements  that  will  be  built  are  initially  iden¬ 
tified  as  air  elements. 

16-20  NA2TYP  Material  type  number  assigned  to  interface 

elements  having  the  properties  of  air.  Usually, 
elements  that  will  be  built  are  initially  iden¬ 
tified  as  air  elements. 
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Column 


Variable 


Explanation 


21-25 

NCTYP 

Structural  material  type,  such  as  concrete  or 
sheetpiling. 

26-30 

NB1TYP 

Backfill  material  type  1.  (Refer  to  section  14b, 
Fill  or  Concrete  Placement.) 

31-35 

NB2TYP 

Backfill  material  type  2.  (Refer  to  section  14b, 
Fill  or  Concrete  Placement.) 

4.  BASE  SEPARATION- -ALPHA  METHOD 

The  following  set  of  variables  control  the  model  of  the  loss  of  contact 
between  a  structure  and  its  foundation  using  the  alpha  method.  Note  that  only 
horizontal  interface  elements  may  be  used  along  the  base  of  the  structure.  If 
this  base  separation  routine  is  not  required,  assign  a  value  of  zero  to  NLOOP 
and  proceed  to  Loading  Information  Card  (section  5). 

(a)  FIRST 

CARD 

FORMAT  (215,  20x,  315) 

Column 

Variable 

Explanation 

1-5 

NLOOP 

Maximum  number  of  iterations  using  alpha  method 
for  each  loading  step. 

6-10 

KEYFRC 

Linear  stress  increment  applied  to  interface 
elements  during  transfer  of  shear  stress  to 
adjacent  element. 

=  0 

31-35 

N0ELF1 

First  horizontal  Interface  element  checked  for 
tension. 

36-40 

N0ELF2 

Last  horizontal  interface  element  checked  for 
tension. 

41-45 

NGROUP 

Number  of  groups  of  interface  elements  for  which 
resultant  forces  are  to  be  computed. 

(b)  SECOND  CARD 

FORMAT  (215) 

Column 

Variable 

Explanation 

1-5 

NELGRP(I,1)  First  interface  element  in  group  I. 

6-10 

MCI PDD/T 

\  JL 

1 2)  Last  liiuCrfsco  aXcmcnt  m  ^rcup  X « 

Repeat  until  I  is  equal  to  NGROUP  sets  of 
interface  elements. 
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5.  LOADING  INFORMATION  CARD 


FORMAT  (613,  2x,  15A4) 


One  card  is  supplied  for  each  loading  step.  One  to  three  loading/ 
construction  modes  can  be  included  in  each  loading  step.  The  loading  or  con¬ 
struction  mode  codes  include: 

KCS(NC.I)  _ DESCRIPTION _ 


1  Excavation  (equivalent  nodal  loads  can  be  applied 
in  equal  increments). 

2  Fill  placement  (subroutine  SUBSTP  cannot  be  used  in 
conjunction  with  the  fill  placement  procedures  of 
subroutine  BUILD) . 

3  Seepage  loading  (equivalent  nodal  loads  can  be 
applied  in  equal  increments) . 

4  Deletion  of  bar  element  (force  in  the  element  can 
be  applied  in  equal  increments) . 

5  Installation  of  bar  element  (prestress  force  can  be 
applied  in  equal  increments) . 

6  Boundary  pressure  loading  (equivalent  nodal  loads 
can  be  applied  in  equal  increments). 

7  Temperature  loading  (the  total  temperature  change 
can  be  applied  in  equal  increments). 

8  Support  displacement  (the  total  displacement  can  be 
applied  in  equal  increments). 

9  Concentrated  nodal  loads  (can  be  applied  in  equal 
increments) . 

10  Element  material  type  change. 

As  indicated  in  the  listing,  input  loads,  displacements,  or  temperature 
changes  can  be  analyzed  in  equal  increments  or  substeps.  Subroutine  SUBSTP 
generates  the  increments,  then  MAIN  analyzes  all  increments  prior  to  analyzing 
the  next  load  step.  With  one  exception,  all  loading/construction  modes  that 
can  be  applied  in  increments,  or  substeps,  can  also  be  applied  in  any  combina¬ 
tion  in  any  loading  step.  The  number  of  substeps,  however,  will  be  the  same 
for  all  loading  or  construction  modes  included  in  the  loading  step.  The 
exception  is  temperature  loading;  if  a  temperature  change  is  specified,  and  a 
given  number  of  substeps  is  specified,  then  only  temperature  loading  can  be 
specified  in  the  loading  step--i.e.,  KCS(N,2)  and  KCS(N,3)  must  be  set  equal 
to  zero  or  left  blank.  If  the  number  of  substeps,  NSBSP,  is  equal  to  zero, 
then  temperature  loading  can  be  included  with  other  loading/construction  modes 
in  a  load  step. 
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Since  th >  same  input  format  is  used  and  similar  opera  ions  are  performed  for 
loading  or  construction  modes  8  and  9,  the  following  ru?es  in  the  usage  of 
these  two  codes  apply.  If  only  concentrated  nodal  loads  are  specified,  use 
mode  9.  If  only  support  displacements  are  specified,  use  mode  8.  If  both 
loads  and  displacements  are  specified,  use  mode  8. 


Column _ Variable _ Explanation _ _ 

1-3  KCS(NC,1)  First  loading/construction  mode  code. 

4-6  KCS(NC,2)  Second  loading/construction  mode  code. 

7-9  KCS(NC,3)  Third  loading/construction  mode  code. 

KCS(NC,1),  KCS(NC,2),  and  KCS(NC,3)  can  be  input  in  any 
numerical  order,  but  the  modes  are  processed  in  ascend¬ 
ing  numerical  order:  If  second  and/or  third  loading 
construction  nodes  are  not  required,  then  KCS(NC,2)  and 
KCS(NC,3)  should  be  set  equal  to  zero  or  left  blank. 


10-12 

NUMIT(NC) 

Number  of  iterations  for  the  loading  step. 
NUMIT(NC)  applies  to  each  substep  if  subste.ps 
are  specified.  A  maximum  of  10  iterations  can 
be  specified.  NUMIT(NC j  «*  0  is  the  same  as 
NUMIT(NG)  -  1. 

13-15 

NUMSS(NG) 

Number  of  substeps.  A  maximum  of  10  substeps  can 
be  specified.  NUMSS(NC)  cannot  be  zero. 

16-18 

M0D(1 ,NC) 

Modulus  calculation  code 

=  1  if  a  loading  modulus  is  to  be  calculated. 
=2  if  an  unload- reload  modulus  is  to  be 
calculated. 


=  0  if  the  computer  is  to  decide  the  type  of 
modulus  to  be  calculated.  In  this  case, 
if  the  most  recently  calculated  maximum 
shear  stress  for  an  element  is  less  than 
all  previous  values  of  maximum  shear 
stress,  an  unload-reload  modulus  is 
calculated.  Otherwise,  a  loading  modulus 
is  calculated. 


Input  is  required  here  only  if  NMOD  “  0.  All  material 
types,  other  than  interface  or  bar  elements,  are  given 
one  of  the  above  codes.  If  NMOD  -  0  and  NC  =  0,  as 
might  be  the  case  for  an  analysis  of  initial  stresses, 
MOD (1,1)  is  set  equal  to  zero,  or  the  computer 
decides . 
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Column 


Variable 


Explanation 


19  IPRT(l.NC)  Print  code 

-  0  print  displacements  and  stresses  for 

final  iteration  only. 

-  1  print  interface  stresses  for  each  itera¬ 

tion,  stresses  for  all  elements,  and  dis¬ 
placements  for  final  iteration  only. 

-  2  print  stresses  for  all  elements  for  each 

iteration,  stresses  and  displacements  for 
final  iteration  only. 

-  3  print  stresses  and  displacements 

for  each  iteration. 

20  IPUN(1,NC)  Disc  storage  code  (tape  9) 

-  0  if  no  disc  storage. 

-  1  if  storage  of  analysis  continua¬ 

tion  data  prior  to  implementation 
of  base  separation  by  alpha 
method. 

=2  if  storage  of  analysis  continua¬ 
tion  data  after  completion  of  base 
separation  by  alpha  method. 

21-80  HEDCS(NC)  Description  of  loading  step. 

6.  MODULUS  SPECIFICATION  CARD  FORMAT  (4012) 

This  card  is  required  only  if  NMOD  -  1  and  NC  -  0.  A  card  is  required  for 
each  loading  step,  1  to  NC.  In  this  option,  values  of  the  modulus  specifica¬ 
tion  code  are  specified  for  each  material  type  (and  thus  each  element,  exclud¬ 
ing  bar  and  interface  elements),  regardless  of  the  change  in  maximum  shear 
strain  that  may  have  occurred. 

Column _ Variable _ Explanation _ 

2,4,6,...  M0D(I,NC)  Modulus  calculation  code  for  each  material 

type  (1  to  NUMSOL)  for  the  first  load  step. 

Sepai ate  cards  are  required  for  each  load  step. 
Columns  not  used  can  be  left  blank. 


7.  MATERIAL  PROPERTY  CARDS 

These  cards  are  used  only  for  two-dimensional  elements;  the  first  and  second 
cards  must  be  supplied  for  each  material  type,  excluding  bar  and  interface 
material  types.  The  first  and  second  cards,  as  a  pair,  are  supplied  in  order 
of  material  type  number  N  =  1  to  NUMSOL.  Information  or  properties  not  re¬ 
quired  for  a  material  type  can  be  left  blank. 

(a)  FIRST  CARD  FORMAT  (7D10.5,  110) 
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Column 

Variable 

Explanation 

1-10 

GUI(N) 

Poisson's  ratio  before  failure. 

11-20 

GUF(N) 

Poisson's  ratio  at  failure  (no  greater  than 
0.49). 

21-30 

GAM(N) 

Total  unit  weight  (always  specified,  regardless 
of  drained  or  undrained  material  behavior) . 

31-40 

FR(N) 

Correlation  factor:  ratio  of  measured 
strength  at  failure  to  ultimate  hyperbolic 
strength. 

41-50 

AO(N) 

Coefficient  of  lateral  earth  pressure  at 
rest,  Kq  ,  pertaining  to  effective  stresses. 

51-60 

PHI(N) 

Friction  angle  in  degrees. 

61-70 

XXP(N) 

Exponent  n  in  the  initial  tangent  and 
unload-reload  modulus  expressions.  Its  value  is 
assumed  to  be  independent  of  loading  mode.  For 
a  linear  elastic  material,  n  must  be  0.  For 
saturated  soils  when  PHI(N)  ■=  0,  n  should 
normally  be  0.0001. 

71-80 

IDRAIN(N) 

Material  behavior  code 
=  0  if  undrained. 

*=  1  if  drained. 

(b)  SECOND 

CARD 

FORMAT  (3D10.3,  D20.5,  3D10.0) 

Column 

Variable 

Explanation 

1-10 

HCOEF(N) 

Coefficient  Km  in  the  initial  tangent 
modulus  expression 

Ei  -  wvv" 

11-20 

ULCOEF(N) 

Coefficient  K  in  the  unload-reload 

modulus  expression 

Eur  "  PaKur^3/Pa)n 


21-30 


COHE(N) 


Undrained  strength  or  cohesion. 


Column 


Variable 


Explanation 


31-50 

E(N) 

Tangent  modulus  at  failure  for  isotropic 
anisotropic  nonlinear  (nonelastic)  materials. 
Young's  modulus  for  elastic  materials. 

51-60 

ALPHA (N) 

Coefficient  of  linear  thermal  expansion  for 
structural  element;  zero  or  blank,  otherwise. 

61-70 

EIMN(N) 

Minimum  initial  tangent  modulus  for 
isotropic  nonlinear  (nonelastic)  materials. 

Zero  or  left  blank  for  elastic  materials. 

71-80 

TENS(N) 

Minimum  allowable  value  of  the  minor 
principal  stress  for  isotropic  nonlinear  (non¬ 
elastic)  materials.  If  tensile,  input  as  a  neg 
ative  value.  Zero  or  blank  for  elastic 
materials . 

8.  INTERFACE  PROPERTY  CARDS 


One  pair  of  cards  are  supplied  for  each  interface  material  type,  N  =  1  to 
NUMJT.  If  no  interface  materials,  no  cards  are  required. 

(a)  FIRST  CARD  FORMAT  (6D10.5,  15) 


Column _ Variable 


1-10 

PHJ(N) 

11-20 

RKS(N) 

21-30 

RKN(N) 

31-40 

COJ(N) 

41-50 

FRJ  (N) 

51-60 

TENSJ(N) 

61-65 

IADJMT(N) 

_ Explanation _ 

Interface  friction  angle  in  degrees. 

Shear  stiffness  before  failure. 

Normal  stiffness  before  failure. 

Interface  cohesion. 

Correlation  coefficent  for  shear  dis¬ 
placement  relationship,  FRJ  =0.0  for  bilinear. 

Tensile  strength  of  interface  material. 

Material  number  of  two-dimensional  soil  element 
next  to  the  interface. 
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(b)  SECOND  CARD 


FORMAT  (2D10.5) 


Column _ Variable 


Explanation 


1-10  STFSMN(N)  Shear  stiffness  after  failure. 

11-20  STFNMN(N)  Normal  stiffness  after  failure. 


9.  NODAL  POINT  CARDS  FORMAT  (110,  4D10.2) 

One  card  is  supplied  for  each  node.  The  numbering  of  nodal  points  must  be 
sequential  and  some  of  the  nodes  can  be  omitted.  Those  nodes  omitted  are 
automatically  generated  by  the  program  at  equal  spacings  between  those  speci¬ 
fied.  The  first  and  last  nodes  must  always  be  specified.  Note  that  DP(N)  and 


PP(N)  are 

automatically  generated  in  equal  increments  for  those  nodes  omitted. 

Column 

Variable 

Explanation 

1-10 

N 

Nodal  point  number. 

11-20 

X(N) 

X  coordinate,  positive  to  right. 

21-30 

Y(N) 

Y  coordinate,  positive  upward. 

31-40 

PP(N) 

Pore  pressure  in  head  of  water;  zero  or 
blank,  if  not  specified.  Pore  pressure  must  not 
be  specified  for  undrained  materials  but  must  be 
specified  for  drained  materials. 

41-50 

DP(N) 

Change  in  pore  pressure  in  head  of  water 

for  soil  elements;  change  in  temperature  for 
linear  elastic  structural  material;  zero  or 
blank,  otherwise. 


10.  BOUNDARY  CONDITION  CARDS 

Cards  1  through  8  are  supplied  as  required  to  specify  restraints  of  boundary 
nodes.  If  there  are  no  nodes  restrained  in  the  mode  specified  for  cards  2,  3, 
or  4,  then  the  card  is  not  required.  For  a  given  card,  specified  nodes  N  =  1 
to  NOY,  NOX,  or  NOXY  must  be  in  sequential  order. 


BIO 


(a) 


FIRST  CARD 


FORMAT  (315) 


Column 

Variable 

ExDlanation 

1-5 

NOY 

Number  of  nodal  points  fixed  against 

Y-movement  only. 

6-10 

NOX 

Number  of  nodal  points  fixed  against 

X-movement  only 

11-15 

NOXY 

Number  of  nodal  points  fixed  against  both 

X-  and  Y-movement. 

(b)  SECOND  CARD 

FORMAT  (1615) 

Column 

Variable 

Explanation 

1-5 

IC(N) 

Nodal  point  number  of  the  first  nodal  point 
fixed  against  Y-movement. 

Additional  nodal  points  fixed  against  Y-movement, 

N  -  2  to  NOY,  are  specified  in  the  next  15  five- 
column  fields  and  on  additional  cards  as  required. 

(c)  THIRD 

CARD 

FORMAT  (1615) 

Column 

Variable 

Explanation 

1-5 

IC(N) 

Nodal  point  number  of  the  first  nodal  point 
fixed  against  X-movement. 

Additional  nodal  points  fixed  against  X-movement, 

N  -  2  to  NOX,  are  specified  in  the  next  15  five- 
column  fields  and  on  additional  cards  as  required. 

(d)  FOURTH 

CARD 

FORMAT  (1615) 

Column 

Variable 

Explanation 

1-5 

IC(N) 

Nodal  point  number  of  the  first  nodal  point 

fixed  against  both  X-  and  Y-movement. 

Additional  nodal  points  fixed  against  both  X-  and 
Y-movement,  N  =  2  to  NOXY,  are  specified  in  the 
next  15  five -column  fields  and  on  additional  cards 
as  required. 
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11 .  ELEMENT  CARD 


FORMAT  (615) 


One  card  is  supplied  for  each  interface  or  two-dimensional  element;  bar  ele¬ 
ments  are  not  included  in  this  series  of  cards.  All  interface  elements  are 
supplied  in  sequential  order  first,  followed  by  two-dimensional  elements,  also 
in  sequential  order.  (If  'Build'  is  used;  place  interface  elements  which  will 
be  built  up  at  the  end  of  joint  list,  i.e.,  start  numbering  these  elements 
with  N  -  NUMJT  -  NJTFIL.)  Thus,  interface  elements  must  be  numbered  from  N  « 

1  to  NUMJT,  and  two-dimensional  elements  from  N  »  (NUMJT  +  1)  to  NUMEL. 

Nodal  point  numbers  must  be  specified  consecutively,  processing  counter¬ 
clockwise  around  the  element.  The  first  and  last  nodal  point  numbers  speci¬ 
fied  for  interface  elements  must  have  the  same  coordinates.  Triangular  two- 
dimensional  elements  having  four  different  nodal  point  numbers  may  not  be 
used;  the  first  and  last  point  numbers  of  a  triangular  element  must  be 
identical. 

Element  numbers  in  a  row  may  be  omitted,  in  which  case  the  omitted  elements 
will  be  generated  by  incrementing  the  element  number  N  and  the  nodal  point 
numbers  I,  J,  K,  and  L  by  one  and  by  assigning  the  same  material  type  number 
as  specified  for  the  last  element.  The  first  and  last  elements  in  the  row 
must  be  specified.  If  no  elements  are  omitted,  the  element  numbering  may  be 
done  in  any  order,  provided  all  interface  elements  are  numbered  first. 

Column _ Variable _ Explanation _ 


1-5 

N 

Element  number. 

6-10 

IL.(N ,  1) 

Number  of  nodal 

point  I. 

11-15 

IL(N , 2) 

Number  of  nodal 

point  J . 

16-20 

IL(N, 3) 

Number  of  nodal 

point  K. 

21-25 

IL(N,4) 

Number  of  nodal 

point  L. 

26-30 

IL(N, 5) 

Material  type  number. 

12.  BAR  ELEMENT  CARD  FORMAT  (415,  2D10.7,  2D10.1,  DIO. 5) 

One  card  is  supplied  for  each  bar  element  initially  in  the  mesh  or,  if  a  con¬ 
tinuation  analysis,  added  in  an  incremental  loading  step  of  the  previous  anal¬ 
ysis.  Note  that  for  a  continuation  analysis  this  card  is  not  automatically 
generated.  Elements  are  numbered  sequentially  from  N  =  1  to  NMBAR. 
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Column 

Variable 

ExDlanation 

1-5 

N 

Element  number. 

6-10 

IB(N , 1) 

Number  of  nodal 

point  I. 

11-15 

IB(N,2) 

Number  of  nodal 

point  J . 

16-20 

IB(N, 3) 

Spring  response 

type  code 

-  1 

if  both  compression  and  tension  of 
bar  allowed. 

-  2 

if  only  compression  allowed. 

»  3 

if  only  tension  allowed. 

21-30 

BAR(N, 1) 

cos  a 

31-40 

BAR(N, 2) 

sin  a 

The  sign  of  angle  a  is  determined  as  shown  in 
Figure  6  of  the  Isotropic  SOILSTRUCT  manual. 

Angle  a  is  measured  counterclockwise  from  a 
line  drawn  in  the  positive  x- direction  at  the 
I  node  to  the  line  of  action  connecting  the 
I  node  to  the  J  node.  (See  comments  that 
follow. ) 

41-50  BAR(N, 3)  Prestress  force  in  the  bar  element.  This  only 

inputs  the  force;  the  force  must  be  applied  as 
loads  at  nodes  I  and  J  using  loading/ 
construction  mode  9  in  a  loading  step. 

51-60  BAR(N,4)  Stiffness  of  bar  element.  This  usually 

computed  as  AE/L,  but  the  mesh  length  (distance 
from  node  I  to  node  J)  need  not,  and  usually 
does  not,  correspond  to  the  actual  length. 

61-70  BAR(N,5)  Displacement  of  bar  element  at  activation. 

This  allows  for  a  specified  degree  of  slack  at 
the  strut  connection;  the  bar  will  deform 
BAR(N,5)  prior  to  its  stiffness  being  activated. 

Bar  elements  can  function  as  either  anchors  or  strut  (spring)  supports.  The 
required  parameters  are  dependent  on  the  type  of  bar  element  specified. 

If  a  strut  support  is  specified,  nodal  point  J  is  a  node  fixed  against  x-  or 
y-movemcnt,  depending  on  the  orientation  of  the  strut  being  modeled.  For 
program  storage  efficiency,  the  number  of  node  J  should  be  as  close  as 
possible  to  the  number  of  node  I.  Nodal  point  I  represents  the  point  of  con¬ 
nection  between  the  wall  and  the  actual  strut.  Nodal  points  I  and  J,  then, 
are  not  necessarily  physically  connected,  since  the  element  stiffness  is  input 
independently.  Nodal  point  J  allows  the  force  at  the  J  node  to  be  carried 
into  the  system  as  a  reaction  at  a  fixed  node.  This  is  consistent  with  the 
typical  mesh  representation  of  one-half  of  a  symmetric  excavation.  The  values 
are  input  to  represent  the  line  of  action  of  the  strut  support  and  do  not  need 
to  correspond  to  the  relative  positions  of  the  I  and  J  nodes. 
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If  an  anchor  is  specified,  nodal  points  I  and  J  physically  represent  two  ends 
of  the  anchor  and  must  be  restrained  appropriately.  The  values  of  sin  a  and 
cos  a  must  correspond  to  the  relative  positions  of  the  I  and  J  nodes  repre¬ 
senting  the  ends  of  the  anchor.  Stiffness  is  computed  as  AE/L,  with  L  being 
the  distance  between  nodes  I  and  J,  and  either  A  or  E  altered  to  give  the 
correct  stiffness.  Stiffness  for  an  anchor  or  a  strut  support  is  input  as 
force  per  length  per  length  of  wall. 

For  either  element  type,  specifying  the  prestress  force  does  not  apply  the 
force.  The  concentrated  force  loading/construction  mode  must  be  used  for  this 
purpose.  Thus,  bar  elements  initially  in  the  mesh  cannot  carry  a  prestress 
force,  since  it  is  not  applied  by  a  gravity  turn-on  analysis. 


13.  CONTINUATION  OR  INITIALIZATION  CARDS 

This  input  is  supplied  only  if  INIT  =0;  it  can  be  supplied  from  external  disc 
storage  or  punched  cards  input  from  a  preceding  analysis.  Input  format  is  the 
same  for  both  input  modes.  This  option  is  provided  so  that  a  required  se¬ 
quence  of  loading  steps  can  be  stopped  at  an  intermediate  step,  then  restarted 
from  that  step  without  redoing  the  complete  analysis.  However,  these  cards 
may  also  be  used  to  specify  values  of  particular  variables  in  an  initial  anal¬ 
ysis  without  using  the  gravity  turn-on  procedure  followed  with  INIT  -  1  or  the 
special  procedure  followed  with  INIT  =  2.  Similarly,  particular  parameter 
values  can  be  changed,  if  the  sequence  of  loading  is  stopped,  prior  to  a 
restart  analysis. 

(a)  FIRST  CARD  (ELEMENT  INFORMATION)  FORMAT  (515) 


Column 

Variable 

Exolanation 

1-5 

NUMEL 

Number  of  elements  in  the  mesh,  excluding 
bar  elements ,  but  including  interface  elements . 

6-10 

NUMJT 

Number  of  interface  elements. 

11-15 

NUMBAR 

Number  of  bar  elements, including  those 
initially  in  the  mesh  and  those  added  in 
previous  loading  steps  (if  a  restart  analysis), 

16-20 

NUMNP 

Number  of  nodes. 
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Column 


Variable 


Explanation 


21-25  NSTART  Interface  element  code 

-0  if  all  interface  elements  are 

active  (i.e.  no  shear  or  tension 
failures).  IFLOLD  and  IFAIL  keys 
are  set  equal  to  zero  for  all 
interface  elements  within  the 
program.  See  card  groups  g  and  h. 
-I  if  IFLOLD  and  IFAIL  keys  are  to  be 
read  for  all  interface  elements. 


When  NSTART  -  1,  the  keys  IFLOLD  and  IFAIL  have 
usually  been  determined  in  a  previous  analysis 
and  are  included  as  standard  output  to  tape  9 
when  the  disc  storage  option  (see  LOADING  INFOR¬ 
MATION  CARD)  was  activated. 

(b)  SECOND  CARD  (STRESS  INFORMATION) 

INTERFACE  ELEMENTS  FORMAT  (2(1P, 1D13 . 6 , 2x) ) 

2 -DIM  ELEMENTS  FORMAT  (4(1P,1D13. 6,2x)) 

Three  cards  are  supplied  for  each  interface  element,  in  numerical  sequence 
N  “  1  to  NUMJT.  For  each  interface  element,  the  pair  of  normal  and  shear 
stresses  at  node  I  [SIGI],  the  center  [SIG] ,  and  node  J  [SIGJ]  are  supplied. 

Column _ Variable _ Explanation _ 

1-13  SIG(N, 1)  Normal  stress  for  an  interface  element. 

16-28  SIG(N,2)  Shear  stress  for  an  interface  element. 

One  card  is  supplied  for  each  two-dimensional  element,  in  numerical  sequence  N 
-  (NUMJT  +  1)  to  NUMEL. 


Column 


Variable 


Explanation 


1-13 


16-28 


SIG(N.l) 


SIG(N, 2) 


Stress  in  the  x-direction  for  a  two- 
dimensional  element,  ax  ,  normal  stress  for  an 
interface  element. 


Stress  in  the  y-direction  for  a  two- 


dimensional  element, 
interface  element. 


shear  stress  for  an 
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Column _ Variable _ Explanation _ 

31-43  SIG(N,3)  x-y  shear  stress  for  a  two-dimensional 

element  r  ,  zero  or  blank  for  interface 
element .  y 

46-58  SIG(N,4)  Maximum  previous  value  of  x-y  shear  stress 

for  a  two-dimensional  element,  zero  or  blank  for 
interface  element. 

(c)  THIRD  CARD  (NODAL  POINT  INFORMATION)  FORMAT  (3(1P1D14 . 7 , lx) ) 

Information  for  three  nodal  points  is  supplied  on  each  card.  Nodal  points  are 
specified  in  numerical  order,  N  =  1  to  NUMNP. 


Column 

Variable 

Explanation 

1-14 

DISPX(N) 

Displacement  in  x-direction 

16-29 

DISPY(N) 

Displacement  in  y-direction 

31-44 

PP(N) 

Pore  pressure  in  head  of  water. 

(d)  FOURTH  CARD  (MATERIAL  TYPE  DESIGNATION)  FORMAT  (1515) 


Materials  type  numbers  for  15  interface  or  two-dimensional  elements  are  speci¬ 
fied  on  each  card.  Material  type  numbers  for  elements  in  numerical  sequence, 

N  *=  1  to  NUMEL,  are  specified. 


Column 

Variable 

Explanation 

1-5 

IL(N,5) 

Material  type  number. 

Material  type  numbers  for  the  next  14  elements 
are  supplied  in  the  next  fourteen  5 -column 
fields . 

Note  that  material  type  numbers  supplied  on  these  cards  supersedes  the  mate¬ 
rial  type  numbers  specified  on  the  element  card  (section  9).  Thus,  material 
type  changes  can  be  made  as  part  of  a  restart  analysis  rather  than  including 
such  changes  in  a  loading  step  of  an  analysis. 

(e)  FIFTH  CARD  (BAR  ELEMENT  INFORMATION)  FORMAT  (4(D10.2,  D10.0)) 

This  card  is  supplied  only  if  bar  elements  are  included  (NUMBAR  >  0) . 


B16 


Information  for  four  bar  elements  is  specified  on  each  card.  Information  is 
supplied  for  bar  elements  in  numerical  sequence,  N  -  1  to  NUMBAR. 

Column _ Variable _ Explanation _ 

1-10  BAR(N,3)  Force  in  bar  element. 

11-20  BAR(N,4)  Stiffness  of  bar  element. 

Information  for  the  next  three  bar  elements  is 
supplied  in  the  next  six  10-column  fields. 

Note  that  these  parameters,  if  changed  for  a  restart  analysis,  supersede  those 
specified  on  the  Bar  Element  Card  (section  12).  Also,  if  bar  elements  are 
initially  included  in  the  mesh,  and  an  initialization  procedure  is  used,  then 
this  card  must  be  included,  duplicating  the  information  specified  in  the  Bar 
Element  Card. 

(f)  SIXTH  CARD  (INTERFACE  INFORMATION)  FORMAT  (8D10.4) 

This  card  is  supplied  only  if  interface  elements  are  included  (NUMJT  >  0) . 
Information  for  four  interface  elements  is  specified  on  each  card.  Informa¬ 
tion  is  supplied  for  interface  elements  in  numerical  sequence,  N  ■  1  to  NUMJT. 
If  no  interface  elements  are  present,  proceed  to  group  12  input  data. 


Column 

Variable 

Explanation 

1-10 

STFS(N) 

Shear  stiffness  of  first 

interface  element. 

11-20 

STFN(N) 

Normal  stiffness  of  first  interface 
element. 

Information  for  the  next 
elements  is  supplied  in 
fields . 

three  interface 
the  next  six  10 -column 

Note  that  the  value  of  the  shear  stiffness,  if  changed  for  a  restart  analysis, 
supersedes  the  value  specified  on  the  Interface  Property  Card  (section  8) . 
Thus,  the  interface  stiffness  can  be  changed  as  part  of  a  restart  analysis. 

(g)  SEVENTH  CARD  (INTERFACE  KEY)  FORMAT  (1515) 

This  card  is  supplied  if  interface  elements  are  included  (NUMJT  >  0)  and 
NSTART  -  1.  Information  for  15  interface  elements  is  specified  on  each  card. 
Information  is  supplied  for  interface  elements  in  numerical  sequence,  N  =  1  to 
NUMJT.  If  no  interface  elements  are  present  or  NSTART  «=  0,  proceed  to  group 
12  input  data. 
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Column 


Variable 


Explanation 


1-5  IFLOLD(N)  Modulus  assignment  key 

-  0  if  interface  element  N  is  active. 

-  1  if  tension  failure  occurred. 

-2  if  shear  failure  occurred. 

(h)  EIGHTH  CARD  (INTERFACE  KEY)  FORMAT  (1515) 

This  card  is  supplied  only  if  NSTART  «=  1  and  for  only  those  interface  elements 
along  the  horizontal  base  of  the  structure  used  to  model  the  base  separation, 
with  element  numbers  from  N0ELF1  to  N0ELF2.  Information  for  15  interface 
elements  is  specified  on  each  card. 

Column _ Variable _ Explanation _ 

1-5  IFAIL(N)  Stress  transfer  key 

“0  if  interface  element  number  N  is 
active . 

*=1  if  transfer  of  normal  and  shear 

stress  in  failed  interface  element 
number  N  has  been  completed. 

14.  LOADING  STEP  CARDS 

These  cards  are  required  only  if  NC  is  not  equal  to  0.  Cards  are  assembled  in 
the  order  specified  on  the  Loading  Information  Card  (section  5):  cards  for 
the  first  loading/construction  mode  specified  for  the  first  loading  step  are 
followed  by  the  cards  for  the  second  loading/construction  mode  specified  for 
the  first  loading  step,  and  so  on  to  the  cards  for  the  last 
loading/construction  mode  specified  for  the  last  construction  step. 

For  a  given  loading  step,  the  lowest  numbered  loading/construction  mode  is 
processed  first,  but  the  analysis  of  the  loading  step  is  made  for  the  combined 
effect  of  all  loading/construction  modes  included  in  the  loading  step.  Care 
must  therefore  be  exercised  in  specifying  some  loading/construction  modes, 
such  as  material  type  changes  or  concentrated  forces  representing  prestress 
forces,  in  the  same  loading  step  with  other  loading/construction  modes. 

(a)  EXCAVATION 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  «*  1. 
Input  is  handled  by  subroutine  EXCAV.  Free  excavated  elements  and 
common  excavated  element  are  input  separately.  A  free  excavated  element 
is  an  element  specified  to  be  excavated  that  has  no  boundary  in  common 
with  an  element  not  specified  to  be  excavated  in  the  loading  step.  A 
common  excavated  element,  therefore,  has  at  least  one  boundary  in  common 
with  an  unexcavated  element. 

Interface  elements  can  only  be  included  as  free  excavated  elements,  even 
if  they  have  a  boundary  in  common  with  an  unexcavated  element.  Inter¬ 
face  elements  cannot  be  used  as  interpolation  elements.  Free  excavated 
elements  (other  than  interpolation  elements)  can  be  used  as  interpola¬ 
tion  elements,  though  common  excavated  elements  are  more  commonly  used. 
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If  possible,  adjacent  common  excavated  elements  should  be  input  sequen¬ 
tially  as  this  procedure  avoids  repetitive  computation;  nodal  loads  for 
a  nodal  point  common  to  the  two  sequential  elements  will  only  be  calcu¬ 
lated  once. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT  (215) 

Column _ Variable _ Explanation _ _ 

1-5  NFXEL  Number  of  free  excavated  elements. 

6-10  NXELCB  Number  of  common  excavated  elements. 

(2)  SECOND  CARD  (FREE  EXCAVATED  ELEMENT  DATA)  FORMAT  (1615) 

Element  numbers  of  16  free  excavated  elements  can  be  supplied  on  one  card.  A 
maximum  of  50  can  be  specified  in  one  loading  step.  Element  numbers  of  all 
free  excavated  elements,  N  «  1  to  NFXEL  are  to  be  specified. 

Column _ Variable _ Explanation _ 

1-5  LNXEL(N)  Element  number  of  first  free  excavated 

element . 

Information  for  the  next  15  free  excavated  ele¬ 
ments  is  supplied  in  the  next  fifteen  5 -column 
fields. 

(3)  THIRD  CARD  (COMMON  EXCAVATED  ELEMENT  DATA)  FORMAT  (815) 

One  card  is  supplied  for  each  common  excavated  element,  N  -  1  to  NXELCB.  A 
maximum  of  50  common  excavated  elements  can  be  specified  in  one  loading  step. 
Loading  codes  include: 

0  -  The  node  is  not  loaded  by  excavation  forces  and  is  not 
common  to  both  an  excavated  and  an  unexcavated  element. 

1  -  The  node  is  loaded  by  excavation  forces  and  is  common  to 
both  an  excavated  and  an  unexcavated  element. 

Note:  I,  J,  K,  and  L  refer  to  the  same  nodes  I,  J,  K,  and  L  specified  on  the 
Element  Card  (section  11) . 

Interpolation  elements  should  be  specified  in  a  crisscross  fashion.  Further, 
x-  and  y-coordinates  of  diagonal  elements  must  not  be  the  same.  If  these 
rules  are  not  adhered  to,  the  interpolation  routine  will  detect  a  singularity 
and  processing  will  stop. 
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Column 

Variable 

Exolanation 

1-5 

LUL(N.l) 

Element  number  of  the  first  common  excavated  ele 
ment.  This  element  is  also  the  first 
interpolation  element. 

6-10 

LUL(N, 2) 

Element  number  of  second  interpolation  element. 

11-15 

LUL(N, 3) 

Element  number  of  third  interpolation  element. 

16-20 

LUL(N,4) 

Element  number  of  fourth  interpolation  element. 

21-25 

LUL(N, 5) 

Loading  code  for  node  I. 

26-30 

LUL(N, 6) 

Loading  code  for  node  J. 

31-35 

LUL(N,7) 

Loading  code  for  node  K. 

36-40 

LUL(N, 8) 

Loading  code  for  node  L. 

Loading  codes  are  specified  for  the  nodes  of  the 
element  to  be  excavated  specified  in  columns  1 
through  5. 

(b)  FILL  OR  CONCRETE  PLACEMENT 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  -  2. 
Input  is  handled  by  subroutine  BUILD.  The  types  of  elements  that  can  be 
placed  include  structural,  soil,  and  interface  elements.  Subroutine 
SUBSTP  cannot  be  used  in  conjunction  with  subroutine  BUILD. 

(1)  F^ST  CARD  (CONTROL  DATA)  FORMAT  (515,  DIO. 2) 


Column 

Variable 

Exolanation 

1-5 

NLEL 

Total  number  of  elements  to  be  placed, 
including  interface  and  structural  elements. 

6-10 

NJ 

Total  number  of  elements  to  be  placed  less  the 
number  of  inactive  interface  elements  to  be 
placed. 

11-15 

NONP 

Number  of  nodal  points  within  the  placed  layers 
to  be  assigned  zero  displacements.  This  in¬ 
cludes  all  nodal  points  of  the  elements  to  be 
placed  except  those  nodal  points  in  common  with 
an  already  existing  element. 

16-20 

NCE 

Number  of  structural  elements  to  be  placed. 
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Column 

Var  .able 

Explanation 

21-25 

NB1E 

Number  of  backfill  type  1  elements  to  be 
placed.  (The  number  of  NB2TYP  elements  is 
determined  internally.) 

26-35 

HTB 

New  y-coordinate  of  the  top  of  the  placed 
backfill . 

(2)  SECOND  CARD  (ELEMENT  NUMBERS)  FORMAT  (1615) 

Element  numbers  of  16  placed  elements  can  be  supplied  on  one  card.  A  maximum 
of  100  can  be  specified  in  one  loading  step.  Element  numbers  of  all  placed 
elements,  N  -  1  to  NLEL,  are  to  be  supplied.  The  elements  must  be  read  in  by 
material  type  according  to  the  following  sequence: 

1.  structural  elements 

2.  backfill  type  1  elements 

3.  backfill  type  2  elements 

4.  interface  elements  to  be  activated. 

5.  interface  elements  to  be  left  inactive  but  to  be  placed 
between  elements  of  like  materials 

Column _ Variable _ Explanation _ 


1-5 

LEL(N, 1) 

Element  number  of  first  element  to  be 
placed. 

6-10 

LEL(N, 2) 

Material  number  of  the  element. 

Information  for  the  next  seven  placed  elements 
is  supplied  in  the  next  fourteen  5 -column 
fields. 

(3)  THIRD  CARD  (NODE  NUMBERS)  FORMAT  (1615) 

Nodal  point  numbers  of  16  nodes  to  be  assigned  zero  displacement  can  be  sup¬ 
plied  on  one  card.  A  maximum  of  50  can  be  specified  in  one  loading  step. 
Nodal  point  numbers,  N  -  1  to  NONP,  are  to  be  specified  in  sequential  order. 

Column _ Variable _ Explanation _ 

1-5  NP(N)  Nodal  point  number  of  first  node  assigned 

zero  displacement. 

Information  for  the  next  15  nodal  points  is  sup¬ 
plied  in  the  next  fifteen  5-column  fields. 


(c)  SEEPAGE 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  «  3. 
Input  is  handled  by  subroutine  SEEP.  Seepage  loads  are  determined  from 
changes  in  pore  pressure  specified  as  DP(N)  on  Nodal  Point  Cards  or  from 
the  specified  phreatic  level  changes . 


B21 


(1)  FIRST  CARD  (CONTROL  DATA) 


FORMAT  (15) 


Column 

Variable 

Exolanation 

1-5 

NCODE 

Option  code  specifying  how  seepage  is 
input 

-0  if  specified  as  DP(N)  on  Nodal 
Point  Cards  (Section  9). 

-  1  if  to  be  calculated  using  the  new 
phreatic  surface  input  on  the  fol¬ 
lowing  cards. 

(2)  SECOND 

CARD  (NUMBER  OF  PHREATIC  SEGMENTS)  FORMAT  (15) 

This  card  is 

required  only  if  NCODE  -  1. 

Column 

Variable 

Explanation 

1-5 

NWAT 

num 

NWAT. 

Number  of  phreatic  surface  segment  end 
points  used  to  specify  the  new  phreatic  surface. 
NWAT  must  be  greater  than  or  equal  to  2.  The 
ber  of  phreatic  surface  segments  is  equal  to 

The  maximum  value  of  NWAT  is  30. 

(3)  THIRD  CARD  (PHREATIC  LEVEL  DATA)  FORMAT  (6D10.2) 

This  card  is  required  only  if  NCODE  -  1. 

The  end  points  of  the  phreatic  surface  segments  delineating  the  new  and  old 
phreatic  surfaces  are  specified  as  x-coordinates  that  must  be  the  same  as  the 
x- coordinates  of  a  nodal  point.  Both  the  present  and  new  phreatic  surfaces 
are  assumed  to  be  linear  between  the  bounding  x-coordinates.  The  left-hand 
side  of  the  mesh  is  always  the  first  bounding  x-coordinate  specified.  A 
bounding  x-coordinate  on  the  old  phreatic  surface  will  require,  usually,  spec 
ification  of  the  same  x-coordinate  on  the  new  phreatic  surface. 

Two  end  points  (x-coordinates),  with  associated  new  and  old  phreatic  levels 
(y-coordiantes) ,  are  supplied  on  each  card.  All  end  points,  N  -  1  to  NWAT, 
must  be  specified. 


Column _ Variable _ Explanation 


1-10 

XW(N) 

X-coordinate  bounding  the  levels  PREL(N) 

and  FUEL(N)  on  the  right-hand  side.  (Must  be 

the  same  as  the  x-coordinate  of  a  nodal  point.) 

11-20 

PREL(N) 

Present  level  (y-coordinate)  of  the  phreatic 
surface  at  XW(N). 

21-30 

FUEL(N) 

New  level  (y-coordinate)  of  the  phreatic 
surface  at  XW(N). 

Information  for  the  next  end  point  is  supplied 
in  the  next  three  10-column  fields. 
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(d)  DELETION  OF  BAR  ELEMENTS 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  -  4. 
Input  is  handled  by  the  main  program. 

The  deleted  bar  elements  remain  in  the  mesh  but  with  zero  stiffness. 
The  force  the  bar  element  carried  is  applied  to  the  free  node  or  nodes 
at  its  ends. 


This  loading/construction  mode  cannot  be  specified  in  the  same  loading 
step  as  fill  or  concrete  placement. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT  (15) 


Column 

Variable 

Explanation 

1-5 

NCARDS 

Number  of  deleted  bar  elements.  There  is 
no  limit  other  than  the  number  of  bar  elements 
presently  in  the  mesh. 

(2)  SECOND  CARD  (ELEMENT  NUMBERS)  FORMAT  (1615) 


The  element  numbers  of  16  deleted  bar  elements  can  be  specified  on  one  card. 
A  total  of  N  -  1  to  NCARDS  cards  must  be  supplied. 


Column 

Variable 

Explanation 

1-5 

N 

Element  number  of  bar  element  to  be 
deleted. 

Element  numbers  for  the  next  15  elements  are 
supplied  in  the  next  fifteen  5-column  fields. 

(e)  ADDITION  OF  BAR  ELEMENTS 


These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  -  5. 
Input  is  handled  by  the  main  program. 

Information  on  the  second  card  is  the  same  as  that  explained  for  the  Bar 
Element  Card  (Section  12) .  The  added  bar  elements  are  numbered  sequen¬ 
tially  from  NUMBAR  +  1,  where  NUMBAR  is  the  number  of  bar  elements  in 
the  mesh  before  the  present  loading  step. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT  (15) 

Column _ Variable _ Explanation _ 

1-5  NCARDS  Number  of  bar  elements  to  be  added.  Any 

number  can  be  added  in  a  loading  step,  however, 
the  maximum  number  of  bar  elements  that  can  be 
in  the  mesh  (including  inactive  or  deleted 
elements)  is  15. 
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(2)  SECOND  CARD 

(ADDED  BAR  ELEMENT  DATA)  FORMAT  (415,  2D10.7,  2D10.1,  DIO. 5) 


Column 

Variable 

Exolanation 

1-5 

N 

Element  number  of  added  bar  element. 

6-10 

IB(N.l) 

Number  of  nodal  point  I. 

11-15 

IB(N,2) 

Number  of  nodal  point  J. 

16-20 

IB(N, 3) 

Spring  response  type  code. 

21-30 

BAR(N, 1) 

cos  a 

31-40 

BAR(N,2) 

sin  a 

Column 

Variable 

Exolanation 

41-50 

BAR(N, 3) 

Prestress  force  in  the  bar  element. 

51-60 

BAR(N,4) 

Stiffness  of  bar  element. 

61-70 

BAR(N, 5) 

Displacement  of  bar  element  at  activation. 

(f)  BOUNDARY  PRESSURE  LOADING 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  -  6. 
Input  is  handled  by  subroutine  3URFLD.  Linear  pressure  distributions 
are  assumed,  based  on  the  pressures  specified  for  the  nodal  points.  The 
right  (positive  x-direction)  and  positive  up  (positive  y-direction) . 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT  (215) 


Column 

Variable 

Exolanation 

1-5 

NLDS 

Number  of  loaded  boundaries.  There  is  no 
limit  to  the  number  of  loaded  boundaries  that 
can  be  specified. 

6-10 

NPRTB 

Print  code 

-  1  if  equivalent  nodal  point  forces 
are  to  be  printed. 
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(2)  SECOND  CARD  (LOADED  BOUNDARY  DATA)  FORMAT  (215,  4D10.2) 


Information  for  one  loaded  boundary  is  specified  on  each  card.  Nodes  I  and  J 
are  interchangeable;  there  is  no  convention  regarding  their  assignment.  A 
total  of  N  -  1  to  NLDS  loaded  boundaries  must  be  specified. 


Column 

Variable 

Evolanation 

1-5 

I 

Nodal  point  number  of  the  first  node  of  the 
load'd  boundary. 

6-10 

J 

Nodal  point  number  of  the  second  node  of 
the  loaded  boundary. 

11-20 

WS1 

Value  of  the  pressure  acting  in  the 
x-direction  at  node  I. 

21-30 

WS2 

Value  of  the  pressure  acting  in  the 
x-direction  at  node  J. 

Column 

Variable 

Exolanation 

31-40 

WS3 

Value  of  the  pressure  acting  in  the 
y-direction  at  node  I. 

41-50 

WS4 

Value  of  the  pressure  acting  in  the 
y-direction  at  node  J. 

(g)  TEMPERATURE  LOADING 

No  cards  are  required  for  this  loading/construction  mode.  Temperature 
Card  (section  9).  If  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  -  7,  then  the 
values  of  DP(N)  are  acknowledged  by  the  main  program  and  processed  as 
temperature  changes.  Note  that  DP(N)  can  also  be  used  to  input  phreatic 
level  changes  for  the  seepage  loading/construction  mode.  Thus,  if  seep¬ 
age  is  specified  as  being  input  through  values  of  DP(N) ,  seepage  and 
temperature  loading  cannot  be  included  in  the  same  loading  step.  Gener¬ 
ally,  temperature  loading  requires  a  restart  analysis,  with  the  DP(N) 
values  being  changed  to  reflect  the  temperature  changes  prior  to  the 
analysis. 

The  temperature  scale  used  (°C  or  °F)  must  correspond  to  the  coefficient 
of  thermal  expansion  designated  on  the  Material  Property  Cards  (sec¬ 
tion  7).  Temperature  changes  are  typically  designated  for  structural 
materials  only. 

(h)  CONCENTRATED  FORCE  OR  DISPLACEMENT  LOADING 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  **  8  or 
9.  Input  is  handled  by  the  main  program.  Refer  to  the  Loading  Informa¬ 
tion  Card  (section  5)  for  instructions  on  using  loading/construction 
modes  8  and  9 . 
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(1)  FIRST  CARD  (CONTROL  DATA) 


FORMAT  (15) 


Column _ Variable _ Explanation _ 

1-5  NUMNDE  Number  of  loaded  or  displaced  nodes.  There 

is  no  limit  to  the  number  of  loaded  or  displaced 
nodes  that  can  be  specified. 

(2)  SECOND  CARD 

DISPLACEMENT  DATA  FORMAT  (2(I5,2(lx,lD14.7))) 

[LOAD  DATA]  FORMAT  (2(I10,2D10.2)) 

Information  for  two  loaded  or  displaced  nodal  points  is  supplied  on  each  card. 
A  total  of  N  -  1  to  NUMNDE  nodes  must  be  specified.  Sign  convention  is  posi¬ 
tive  to  the  right  (positive  x-direction)  and  positive  up  (positive 
y-direction) .  Nodal  points  specified  as  being  loaded  or  displaced  do  not  have 
to  be  in  numerical  order. 


Column _ Variable 


Explanation 


1-5 

[1-10] 

I 

7-20 

[11-20] 

XI 

22-35 

[21-30] 

Y1 

36-40 

[31-40] 

J 

42-55 

[41-50] 

X2 

57-70 

[51-60] 

Y2 

Node  number  of  the  first  loaded  or 
displaced  node. 

Component  of  force  or  displacement  in  the 
x-direction  at  node  I. 

Component  of  force  or  displacement  in  the 
y-direction  at  node  I. 

Node  number  of  the  second  loaded  or 
displaced  node. 

Component  of  force  or  displacement  in  the 
x-direction  at  node  J. 

Component  of  force  or  displacement  in  the 
y-direction  at  node  J. 


If  there  is  not  a  second,  or  J,  node  to  be  specified  on  the 
last  card,  then  leave  the  second  set  of  columns  blank. 


.'i)  ELEMENT  MATERIAL  TYPE  CHANGE 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  =  10. 
Input  is  handled  by  the  main  program. 

The  material  type  of  the  specified  element  is  changed  before  the  analy¬ 
sis  of  the  loading  step  which  specifies  the  change.  The  material  type 
change  includes  modifying  the  values  of  modulus,  E,  and  Poisson's  ratio, 
GUI,  and  zeroing  the  stresses,  SIG(N,1).  Thus,  if  a  material  type 
change  is  specified  in  conjunction  with  boundary  loading  in  the  same 
loading  step,  the  elements  with  changed  material  type  will  respond  to 
the  loading  with  new  material  properties. 
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As  included,  this  loading/construction  mode  is  intended  to  physically 
represent  the  grouting  of  an  anchor.  At  a  given  step  in  the  analysis, 
the  material  types  of  soil  elements  can  be  changed  to  represent  the  as¬ 
sumed  linear  elastic  grout  zone.  If  there  is  a  need  to  change  material 
types  for  any  other  reason,  this  can  be  done  by  stopping  the  analysis 
after  the  appropriate  lo?d  step,  modifying  the  material  types  on  the 
Material  Type  Designation  Card  (section  13d) ,  then  restarting  the 
analysis . 

(1)  FIRST  CARD  (CONTROL  CARD)  FORMAT  (15) 

The  maximum  number  of  elements  with  material  type  number  that  can  be  changed 
in  a  load  step  is  120.  The  excavation  and  material  type  change 
loading/construction  modes  cannot  be  specified  in  the  same  loading  step  since 
the  same  variable,  LUL(N,1)  and  LUL(N,2),  are  used  to  input  data  for  both. 

Column _ Variable _ Explanation _ 

1-5  NELCH  Number  of  elements  with  material  type 

number  being  changed. 

(2)  SECOND  CARD  (ELEMENT  DATA)  FORMAT  (1615) 

The  element  numbers  and  new  material  type  numbers  of  eight  elements  can  be 
supplied  on  one  card.  A  total  of  N  =  1  to  NELCH  elements  must  be  specified. 

Column _ Variable _ Explanation _ 

1-5  LUL(N,1)  Element  number  of  first  element  with  a 

specified  new  material  type  number. 

6-10  LUL(N,2)  New  material  type  number  of  the  specified 

element. 

Information  for  the  next  seven  elements  is  supplied  in  the 
next  seven  pairs  of  five-column  fields. 
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